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ABSTRACT 

Following our study on the incidence, morphology and kinematics of the ionised gas in 
early-type galaxies we now address the question of what is powering the observed nebular 
emission. To constrain the likely sources of gas excitation, we resort to a variety of ancillary 
data, we draw from complementary information on the gas kinematics, stellar populations 
and galactic potential from the SAURON data, and use the SAURON-specific diagnostic dia- 
gram juxtaposing the [O 111] A5007/H/3 and [N l] AA5197, 5200/H/3 line ratios. We find a tight 
correlation between the stellar surface brightness and the flux of the H/? recombination line 
across our sample, which points to a diffuse and old stellar source as the main contributor of 
ionising photons in early-type galaxies, with post-asymptotic giant branch (pAGB) stars being 
still the best candidate based on ionising-balance arguments. The role of AGN photoionisa- 
tion is confined to the central 2" — 3" of an handful of objects with radio or X-ray cores. 
OB-stars are the dominant source of photoionisation in 10% of the SAURON sample, whereas 
for another 10% the intense and highly-ionised emission is powered by the pAGB population 
associated to a recently formed stellar subcomponent. Fast shocks are not an important source 
of ionisation for the diffuse nebular emission of early-type galaxies since the required shock 
velocities can hardly be attained in the potential of our sample galaxies. Finally, in the most 
massive and slowly- or non-rotating galaxies in our sample, which can retain a massive X-ray 
halo, the finding of a spatial correlation between the hot and warm phases of the interstellar 
medium suggests that the interaction with the hot interstellar medium provides an additional 
source of ionisation besides old UV-bright stars. This is also supported by a distinct pattern 
towards lower values of the [O lll]/H/3 ratio. These results lead us to investigate the relative 
role of stellar and AGN photoionisation in explaining the ionised-gas emission observed in 
early-type galaxies by the Sloan Digital Sky Survey (SDSS). By simulating how our sample 
galaxies would appear if placed at further distance and targeted by the SDSS, we conclude 
that only in very few, if any, of the SDSS galaxies which display modest values for the equiva- 
lent width of the [O ill] line (less than --2.4A) and LINER-like [O lll]/H/3 values, the nebulai- 
emission is truly powered by an AGN. 
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1 INTRODUCTION 



If more than 30 years of investigations have now asserted the 
common presen ce of ionised gas in early-type galaxies (e.g . 



Caldwell 1984; Phillips et al. 1986; Kim 1989; Shields 1991; 
Goudfrooii et al, 1994. ; Macchetto et al. 1996 ; Sarzi et al. 2006) , 



the sources powering the observed nebular emission have still to 
be firm ly identified. As alrea dy revealed by the spectroscopic sur- 
vey of IPhillips etai] ( Il986h . in early-type galaxies the ionised- 
gas emission show values for several line ratios that are gener- 
ally consiste nt to what obser ved in low-ionisation nuclear regions, 
or LINERs ( lHeckman|[l98(]h . Since the na ture o f LINER activity 
in itself is somehow controversial (see^Hol llOOSl , for a review), it 
is perhaps no surprise that the excitation mechanism for the ex- 
tended LINER-like emission is still not known. In fact, in addi- 
tion to the two most touted sources for LINER activity, a central 

AGN or fast shocks, the interaction with the hot phase of the inter- 

J I I 

stellar medium dSparks. Macchetto, & Golomb ek 1989) and pho- 

toionisation by old UV-bright stellar sources (Binette et al. 1994) 
have also been suggested as mechanisms for powering the nebu- 
lar emission of early-type galaxies. The importance of these last 
two ionising source wa s highlighted by the narrow-band survey of 
iMacchetto et al. (1996) and further stressed by the spectroscopic 
observations of Goudfrooij ( 199^) , albeit only in the case of giant 
and dusty elliptical galaxies. 

In this context , the SAURON integral-field survey 
dde Zeeuw et"ai] l2002h allowed a more systematic spectro- 
scopic study of the ionised gas in ne arby early- t ype g alaxies, to 
an unprecedented sensitivity limit. In ISarzi et alj j2006l. hereafter 
Paper V) we reported an amazing variety for the value of the 
[O IIl]A5007/H/3 ratio both across the SAURON sample and within 
single galaxies, which suggests that a number of mechanisms is 
responsible for the gas excitation in early-type galaxies. If more 
than one source contributes to power the extended emission of 
early-type galaxies, disentangling their relative importance in 
nearby objects could prove critical also to understand the gas 
emission observed in distant galaxies by large-scale surveys 
such as the Sloan Digital Sky Survey (SDSS), where the nebular 
fluxes are integrated over large physical apertu res. Yet, the limited 
wavelength range of the SAURON data (see iBacon et al] l200ll) 
prevents the standard emission-line diagnostic that is routinely 
used to make a first assessment of the poss ible ionising sources 
(e.g., following I Veilleux & Osterbrock|[l987l) . so that other pieces 
of information are needed to tackle this problem. Presently, we can 
add to our emission-line measurements not only our knowledge of 
the g ravitational potential of our sample galaxies ( Capp ellari et al.l 
120061 . Paper IV) and of the basi c properties of their stellar 
populations jKuntschner et alj|2006l . Paper VI), but we can also 
draw from a variety of ancillary data that have been collected and 
analysed over the course of th e SAURON project . These range from 
GALEX ultraviolet images jjeong et alj |2009|) , Spitzer infrared 
data jShapiro etai] 12003) to mm- and cm- radio observations of 
the molecular and n eutral gas dCombes. Young. &Bureau|[2007l ; 
iMorganti et ai] l200^). Furthermore, for a minority of our sample 
galaxies it is possible to use a SAURON-specific diagnostic diagram 
that uses the [N l]/H/3 line ratio as a gauge for the hardness of the 
ionising continuum. We are therefore in a unique position to start 
assessing the relative role of the various sources of ionisation that 
could be responsible to the nebular emission observed in early-type 
galaxies. 

Our investigation will proceed as follows. In ^we will in- 
troduce the potential of the [N l]/H/3 diagnostic diagram and derive 



more robust [N I] measurements in order to maximise its use across 
our sample. In ^we will review the importance of AGN activity, 
star formation, diffuse old stellar sources, fast shocks and of the 
interaction with the hot phase of the interstellar medium in power- 
ing the nebular emission observed across the SAURON sample. In 
Sj4]we will explore the implications of the lessons learnt in ^for 
interpreting the ionised-gas emission observed in distant early-type 
galaxies that are targeted by large-scale spectroscopic survey such 
as the SDSS, and finally draw our conclusions in ^ 



2 SAURON EMISSION-LINE DIAGNOSTICS 

Emission-line diagnostic diagrams have had a long history of suc- 
cess in constraining the different sources of ionisation for the nebu- 
lar emi^sion_observedj^ Initially introduced 
by [Baldwin. Phillips. & TerlevichI ( ll98lL hereafter BPT) to sepa- 
rate star-forming objects and active galactic nuclei (AGN), these 
diagrams juxtapose various emission-line ratios to provide a ba- 
sis for classifying the gas emission into categories that, to vary- 
ing degrees of confidence, can be identified with specific excita- 
tion sources. This is achieved by comparing the position of the 
data in the diagnostic diagrams with the predictions of sophisti- 
cated models for gas that is pho toionised by a central AGN or by 
OB-stars (e.g.. lFerland et aljl998 ; Kewlev et al. 2001) or that is ex- 
cited by shocks (e.g.. lDopita & Su therland 1995i , i 19961 . Today, the 
most widely a dopted kind of BP T-diagrams are those introduced 
bv lVeilleux & Osterbrockl h987h . who designed a diagnostic anal- 
ysis that is insensitive to reddening and that is based on the ratios of 
strong emission lines rather close in wavelength, such as H/3 A4861 
and [O III] AA4959, 5007 or [O I] A6300, [N ii] AA6548, 6583, 
Ha A6563, and [S II] AA6716, 6731. 

Among the emission lines that usually feature in these dia- 
grams only the [O III] and H/3 lines fall in the wavelength range of 
the SAURON spectra, which was intentionally restricted around the 
H/9 and Mg6 absorption-line features in order to measure the stel- 
lar kinematics and the properties of the stellar populations within 
the largest possible field of view. A traditional diagnostic anal- 
ysis in our sample galaxies is therefore only possible when the 
SAURON data are complemented by integral-field observati ons in 
the Hq+[N II] spectral region, as shown bv lMazzuca et alj j2006l) 
for the Sa galaxy NGC 7742. 



2.1 The [N I] lines and the [N i]/H/3 diagnostic diagram 

If standard BPT-diagrams catmot be drawn with SAURON data 
alone, the emission from the weak [N I] AA5197, 5200 doublet al- 
lows for a diagnostic analysis that is rather similar to that normally 
carried out through the BPT-diagram based on the [O I]/Ha line ra- 
tio. Indeed, similarly to the [O I] emission also the [N I] lines arise 
in partially ionised regions, which are extended in gaseous nebu- 
lae photoionised by a spectrum containing a large fraction of high- 
energy photons, but are almost absent in H Il-regions. The [N l]/H/3 
ratio can therefore be used to gauge the hardness of the ionising 
continuum, making the [N l]/H/3 vs. [O IIl]/H/3 diagnostic diagram 
in principle as useful as the [O I]/Ha vs. [O IIl]/H/3 diagram to sep- 
arate photoionisation by OB-stars from other excitation sources. 

This is demonstrated by Fig. [T] which shows that the predic- 
tions of the MAPPING-III models in the event of a starburst, of gas 
photoionised by a central AGN and of gas excited by shocks (from 



Dopita et al .l2000l ; lGroves. Dopita. & Sutherlanj2004l ; l Allen et al.1 
20081 respectivelv) are as well separated in the SAURON [N l]/H/9 
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Figure 1. The standard [O I]/Ho vs. [O III]/H/3 BPT diagram of Veilleux & Osterbrock (left) compared to the SAURON [N I]/H/3 vs. [O in]/H/3 diagnostic 
(right). In both panels the red, blue and green lines show the predictions of the MAPPINGS-III models for gas that is photoionised by a central AGN, by 
O-star s, or that is excited by shocks, respectively. The AGN grids are from the dusty, radiation pressure-dominated models of lGroves. Dopita. & Sutherland 
and adopt three values for the index a of the power-law AGN continuum oc i/" {a = —1.7, —1.4, —1.2 from left to right). In each AGN 
model grid the solid lines trace the dimensionless ionisation parameter log U (defined as log q/c), which increases with the [O III]/H/3 ratio from log U = 
—3.0, —2.6, —2.3, —2.0, —1.6, —1.3, —1.0, whereas the dotted lines show the adopted values for the electron density of Ne = 10^ and 10* cm~^, with 
smaller values of the [N I]/H/3 ratio corresponding to larger Ne values. The s tarburst grids are fromlPopita et al.l i2000l) . and assume a gas density of A^e = 
350 cm~^ and use a spectral energy distribution obtained from Starburst99 jLeitherer et alj[l999l) models for an instantaneous star-formation episode. The 
grids assume a range of metallicity Z for both stars and gas in the starburst, shown by the solid lines for Z = 0.2, 0.4, 1.0, 2.0 Zq, and different values of the 
ionising parameter q, shown by the dotted lines for q = 0.5, 1, 2, 4, 8, 15, 30 X 10^ cm s~^. Similarly, for the shock grids (without precursor H 11 region), the 
solid lines show models with increasing shock velocity Vs = 150 , 200, 300, 500, 750, 1000 km s~^, and the dotted lines models with magnetic parameter 
b = 0.5, 1.0, 2.0, 4.0. The shock giids are from lAllen erai] (2008) and assume an electron density of Ne = 1 cm and a solar value for the gas metalhcity. 
In both panels the filled circles, the crosses and the open circles show th e v alues of the emission -line ratios observed in the central 3" of 284 galaxies with 
SDSS spectra, where according to the limits draw bv lKewlev et alj feoOlh and lKewlev et alj i200d) in the [O I]/Ha diagram the nebular emission is most likely 
arising from H II regions, from a Seyfert AGN or from LINER-like emission regions, respectively. 



diagnostic as they are in the standard [O I]/Ha BPT diagram. In 
Fig. [T] we also show the distribution of the emission-line ratios 
measured with the method of Paper V in the central regions of 
284 galaxies observed by the Sloan Digital Sky Survey (SDSS). 
This sample includes all galaxies with SDSS spectra (from the 
DR6 release), not only early-type galaxies, with redshifts between 
0.02 < z < 0.05, apparent r-band magnitude brighter than r < 16 
and, most important, with detected [N I] emission. If we adopt 
the conservative theoretical limit of iKewlev et alj ( 1200 ll) to iden- 
tify star-bursting SDSS galaxies in the [O \]fHa BPT diagram and 
further separate Seyfert nuclei from objects showin g LINER-like 
centr al emission using the empirical line drawn bv iKewlev et al.l 
( I200 6). Fig.[T]illustrates how these three kinds of objects are clearly 
separated also in the SAURON [N l]/H/9 diagnostic diagram. 

The position of the SDSS galaxies displaying nuclear star for- 
mation and LINER-like emission in Fig. [T] deserves further at- 
tention. Starting with the star-forming galaxies, the fact that the 
MAPPINGS-III starburst models systematically underpredict the 
values of the [N l]/H/3 ratio suggests that, similarly to case of the 
[O I] line ( Stasiiiska & Leitherer 1996; Dopita et al. 2000), also the 
[Nl] fluxes can be significantly enhanced by the mechanical en- 
ergy released by supernovae and stellar winds, which are naturally 
expected in starbursts as a result of stellar evolution. As regards the 
location of the SDSS galaxies exhibiting central LINER-like emis- 



sion, we note that these objects are more clearly separated from 
the Seyfert nuclei in the [N l]/H/3 diagnostic diagram than they are 
in the [O I]/Hq diagram. Furthermore, most of the SDSS LINER 
nuclei appear to follow a tight sequence consistent with photoioni- 
sation by a central AGN with increasingly harder ionising radiation 
fields, in agreement with the findings of Kewley et al. (2006) based 
on the [O I]/Hq diagnostic. These results further stress the useful- 
ness of the [N I] lines as a gauge for the hardness of the ionising 
continuum. 



Unfortunately, the 2"^^ and 3'^'^ energy levels of neutral Nitro- 
gen have a much smaller critical de-excitation density (7 x 10^ and 
2.2 X 10^ cm"^) than the 4* level of neutral Oxygen (1.5 x lO** 
cm~'^), so that the [N I] lines are generally much fainter than the 
[O I] line owing to collisional de-excitation. The SAURON [N l]/H/3 
diagnostic diagram is therefore applicable only in the few early- 
type galaxies where nebular emission is particularly intense. To 
complicate matters further, as discussed in Paper V the detection 
of the [N I] lines is hampered by our limited ability to match the 
underlying continuum in the spectral region around the Mgfe ab- 
sorption features. 
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2.2 Improving the [N I] measurements with the MILES 
spectral library 



In ^Capp ellari et al ] hoot Paper X) we re-extracted the stellar 
kinematics for our sample galaxies using for each galaxy a dif- 
ferent subset of stellar temp lates from the MILES library of 
ISanchez-Blazquez et alj ( l2006h . The MILES templates provide a 
much better fit to the emission-free spectral regions of our data, 
and thus their use delivers a more reliable stellar kinematics. 

In an attempt to further improve the extent and quality of 
our [N I] measurements, we re-extracted also the nebular emission 
fluxes and kinematics using the MILES stellar library. In the case 
of emission-line measurements, however, it is crucial to provide an 
underlying fit to the stellar continuum that is as physically moti- 
vated as possible in order to avoid spurious results. For this reason, 
rather than using a subset of stellar spectra from the MILES library 
for each galaxy as done in Paper X, we used a t emplate library 



Vazdekis et al 
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consisting of simple stellar population models from 
(120091 based on MILES stellar spectra, see also 
uOOTD to which we added a number of empirical templates ob- 
tained by matching real SAURON spectra, most of the times devoid 
of emission. 

Specifically, based on the maps of Paper V for the ionised- 
gas emission, we extracted a number of high-S/N SAURON spectra 
from circular apertures in galactic regions either without significant 
emission or with emission lines that were clearly distinguishable 
from the underlying absorption spectrum, which we fitted over the 
entire wavelength range with the MILES stars (excluding stellar 
templates with exceedingly low values of the H/3 absorption-line 
strength). The optimal combination of the stellar spectra that best 
matched each of our emission-free aperture spectra constitute each 
of the empirical templates in our new template library. Care was 
taken to extract SAURON spectra covering the distribution of the 
absorption-line strengths observed in Paper VI, which cannot be 
covered in detail by the distribution of line-strengths values mea- 
sured in the MILES stellar library since this consists mostly of 
spectra for nearby stars. In this respect, our empirical templates 
approximate as closely as possible the spectra of real early-type 
galaxies, were they unaffected by kinematic broadening. Fig. |2]il- 
lustrates the poor coverage that the MILES stellar library provides 
for the absorption-line strengths observed in early-type galaxies, 
in particular in the plane defined by the Mgfe and Fe5015 indices, 
which the adopted array of single stellar population models and 
empirical templates can match in a more physical manner. 

Compared to the fit of Paper V, the use of such a mixed library 
of model and empirical templates based on MILES stellar spectra 
leads to a significant increase in the quality of the fit to the SAURON 
spectra in the Mg6 region (corresponding to a 30% decrease in the 
RMS of the fit residuals), which adds confidence in the [N I] mea- 
surements. In general the H/3 and [O III] measurements obtained 
with the new set of templates do not differ significantly from our 
previous measurements, which reflect the modest 10% decrease in 
the RMS of the fit residual when the entire SAURON wavelength 
range is considered. Therefore, for the sake of clarity hereafter we 
will quote and plot the H/3 and [O III] measurements from Paper V, 
and use the newly extracted emission-line measurements only when 
discussing the [N l]/H/3 diagnostic diagrams or when comparing 
maps for the [N I] and H/3 emission. 




Figure 2. Distribution of the values for the Lick Mgfe and Fe5015 indices 
observed across the SAURON sample (from the maps of Paper VI, black 
dots) compai'ed to that observed across the MILES stellai' librai'y (left, blue 
circles) or for the single-age models of Vazdekis et al. and the empirical 
templates based on matching real SAURON spectra (right, red and orange 
circles, respectively). The MILES stars provide a poor coverage of the re- 
gion in the Mgfe vs. Fe5015 plane that is most populated by the SAURON 
data, which causes the continuum fit process to adopt unphysical combi- 
nations of stars and thus leads to spurious emission-line measurements, in 
particulai' to exceedingly large [N I] fluxes. This is not the case when simple 
stellar population models and empirical templates are adopted, although the 
quality of the fit is formally not as good. 



3 lONISATION SOURCES 

Except for very low and high values of the [O IIl]/H/3 ratio that 
most likely correspond to H Il-regions and Seyfert nuclei, respec- 
tively, the diagnostic diagrams of Fig.[T]illustrate vividly how with- 
out other line ratios than [O IIl]/H/3 it is impossible to immedi- 
ately separate star formation from AGN activity. Besides, a number 
of ionis ing mechanisms in addition to AGN activity, such as fast 
shocks iDopita & Sutherland! 1 19951) . photoionisation by old stars 
teinette et al.'l994') and the interaction with the hot. X-ray emitting 
gas (Sparks, Macchetto, & Golombek 1989), have been suggested 
to explain the emission in the log([0 IIl] /H/3)= - 0.5 range, which 
is typ ical of LINER-like emission (e.g.. lHo. Filippenko. & SargentI 
Il997h and that includes 75% of [O IIl]/H/3 values observed across 
the SAURON sample. In fact, very low [O IIl]/H/3 values are ob- 
served also in the absence of star formation, for instance in the 
X-ray filaments a r ound the brightest galaxies of cl usters (e.g., 
iHatch et al. I l2006l : ISabra. Shields. & Filippeiikol I2OOQ) , whereas 
quite large [O IIlJ/H/3 ratios can be observed also in H Il-regions 
if the metallicity is extremely low. 

For this reason in this section we will use a variety of ancil- 
lary data, ranging from radio to X-ray wavelengths, and consider 
complementary information on the gas morphology and kinemat- 
ics from our SAURON data to constrain the most likely sources of 
gas excitation in different kind of galaxies, using, when possible, 
the [N l]/H/3 diagnostic to corroborate our findings. 



3.1 Nuclear Activity 

Many of our sample galaxies host AGNs. Well-known cases in- 
clude NGC 4486 and NGC 4374 for the ir FR I radio jets (e.g., 
lOwen et alj 1 19891 : iLaing & Bridie 1 1987b NGC 4278 for host- 
ing a well-studied LI NER nucleus (e.g., iGoudfrooii eTai]|l994l: 
iGiroletti. Taylor. & G iovannini 2005) and NGC 4552 for an ultra- 
violet flare that could have resulted from th e tidal stripping of a 
star passing close to the central black hole dRenzini etalJll995l: 
ICappellariet"ai]| 19991) 

Given that unresolved radio-continuum cores are excellent 
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Figure 3. SAURON H/3 emission vs. VLA measurements for the 1.4 Ghz radio continuum from the FIRST survey and for galaxies with compact radio cores. 
As in Paper V the H/3 fluxes are colour-coded in a logarithmic scale with values ranging between 10"^'* and 10"^** erg s~^cm~^arcsec^^. The red contours 
show the 1.4 Ghz fluxes also in logarithmic scale, starting from a level of 0.45 mjy/beam (corresponding to 3 times the typical RMS of the FIRST images) 
and doubling with each contour, except for NGC 4278, NGC 4374, NGC 4486 and NGC 4552 for which the 1.4 Ghz fluxes changes by a factor 5 between 
each line. The beam of the FIRST images is circular and has a FWHM of 5'.'4. 



signposts for black-hole accretion, in this section we focus on the 
galaxies in our sample that exhib it such a feature. In the area cov- 
ered by the VLA FIRST survey jSecker. White. & Helfan3l 19951) 
there are 15 galaxies with 1.4 Ghz compact cores. Figure [3] com- 
pares the SAURON H/9 emission with the distribution of the 1 .4 Ghz 
radio continuum, which in several cases extends beyond the nu- 
cleus. In particular, a comparison with the beam of the FIRST im- 
ages reveals that the nuclear emission is resolved in NGC 2974, 
NGC 3414, NGC 4546, NGC 5813 and NGC 5846, while in 
NGC 3489 the radio core is only marginally extended. 

Among the SAURON galaxies that have not been observed by 
the FIRST survey, the presence of a radio core has been excluded 
in NGC 474 and N GC 821 bylWrob el ( 1991) and in NGC 7457 
and NGC 7332 by iMorganti etlj] I2OO6I The remaining galax- 
ies, NGC 524 and NGC 1023, do host a radio core according to 
iFilho et alJ ( l2004h and Morganti et al., respectively. In what fol- 
lows we will use the radial profiles for the values of the H/3 flux 
and the [O IIl]/H/3 ratio to understand the extent to which AGN can 
power the nebular emission in our sample galaxies. In particular, 
it will be shown that extremely shallow H/3 and [O IIl]/H/3 are in- 
consistent with AGN photoionisation. This is the case of NGC 524, 
whereas for NGC 1023 the absence of central H/3 emission pre- 
cludes any further analysis. Conversely, among the objects with un- 
detected radio-continumm cores in the FIRST images, NGC 4477 
display a strong and unresolved central H/3 peak and a very steep 
[O IIl]/H/3 gradient, suggesting that AGN photoionisation is impor- 
tant toward the centre of this galaxy. The presence of a powerful 



AGN in NGC 4477 is signalled by a strong X-ray nucleus (first 
detec t by the Einstein observatory bv lFabbiano, Kim. & Trinchieril 
and its radio-quiet nature is consistent with its Seyfert clas- 
sification (Ho et al. 1997:' HoI|2008|) . 

If we assume a uniform gas density and constant filling fac- 
tor for gas clouds residing on a plane, in ionisation equilibrium the 
surface brightness of the H/3 recombination line should follow the 
strength of the radiation field, therefore radially decreasing as 
in the case of a central source. The strength with which the AGN 
radiation reaches the gas clouds, usually quantified by the ionisa- 
tion parameter q that compares the ionising-photon and electron 
density at the face of the irradiated cloud, also determines the ioni- 
sation state of the gas clouds. For our purposes the [O IIl]/H/3 ratio 
is an excellent diagnostic for the ionisation state of the gas, not only 
since this ratio involves emission from a highly-ionised spe cies 
(e.g.. lFerland & Netzerlll983l : lHo. Shields. & FiliDpenkdll993h but 
also because the [O IIl]/H/3 ratio scales almost linearly with q 
for the range of values observed in the SAURON nuclei (i.e., for 
[O IIl]/H/3^ 5). Therefore, also the values of the [O IIl]/H/3 ratio 
should follow a r^^ profile under the previous assumptions. 

On the other hand, if the gas density decreases with radius, 
the H/3 surface brightness will drop at an even faster gradient than 
inverse-square, whereas the [O IIl]/H/3 ratio will display a shal- 
lower radial profile as the ionisation parameter q decreases more 
gently with radius. Furthermore, if the gas is also vertically ex- 
tended rather than just being distributed in a planar configuration, 
projection effects will lead to shallower H/3 surface brightness pro- 
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Figure 4. Radial profiles for the H/3 flux values in galaxies wifli 1.4 Ghz radio cores. The H/3 fluxes are in units of 10~^^erg s~^cm~^arcsec^^ and are 
plotted only out to 15". In each panel the solid hne shows the radial decrement that is expected in the case AGN photoionisation and under the special 
conditions of constant gas density and filling factor, accounting for the impact of atmospheric blurring. The point-spread function of the S AURON observations, 
from Paper 111, is shown by the dotted lines. Both profiles have been normalised to the central H/3 emission. The H/3 data are colour-coded according to the 
logarithm of the [O III]/H/3 ratio, going from blue to green and red for values of -1, and 1. 



files. In this case the [O IIl]/H/3 values will follow the same trend as 
the H/3 emission and also show a gentler radial decline, as clouds 
at different distances from the centre and hence subject to different 
ionising photon densities are observed along the same line of sight. 
Conversely, if the gas clouds lie in a plane, the AGN radiation could 
be absorbed by intervening material before reaching the gas clouds 
at any given radius, leading to H/3 and [O ni]/H/3 profiles that are 
steeper than an r^^ law. Finally, radial gradients in the filling factor 
of the gas clouds will only affect the H/3 surface-brightness profiles, 
in the same direction as variations in the gas density would do. 

Figure |4] plots the radial profile of the H/3 emission in the 
S AURON galaxies with 1.4 Ghz radio cores, together with the pre- 
dicted radial decrement for AGN photo-ionisation and under 
the special conditions of constant gas density and filling factor, ac- 
counting for atmospheric blurring. The H/3 data-points in Fig. |4] 
are colour-coded according to the value of the [O IIl]/H/3 ratio, in 
order to convey also the radial variation of this quantity. Figure |4] 
shows that the surface brightness profile of the H/3 emission is in 
keeping with our simplest expectations for AGN photo-ionisation 
only in the very central regions of a few objects (e.g., NGC 3414, 
NGC 4278, NGC 4546) whereas in general the H/3 profile is shal- 
lower than an inverse-square law, in particular beyond the central 
~ 3". The behaviour of the [O IIl]/H/3 ratio, which does not appear 
to systematically decrease with radius, further rules out a central 
AGN as the source of ionisation for the gas emission beyond the in- 
nermost regions of our sample galaxies. In fact, that the [O IIl]/H/3 



ratio tends to display rather constant values in the innermost re- 
gions (within 20% of the average within 3", see also Fig. 4b of 
Paper V) suggests that if an AGN is responsible for the gas exci- 
tation where the H/3 surface-brightness profile is consistent with a 
r^^ law, then the assumptions of our simplest scenario cannot hold. 
In this respect, we note that if the central gas of our sample galax- 
ies has settled on a steady disk, this should display a radially de- 
creasing profile for the column density ( Kawata, Cen, & Ho 200"^, 
which is proportional to the gas density and the filling factor of the 
gas clouds. The presence of a negative gradient for the gas density 
could address the need for shallow [O IIl]/H/3 profiles in the nuclear 
regions of galaxies with H/3 gradients that are consistent with the 
inverse-squared model, even though a radially decreasing gas den- 
sity would induce even steeper H/3 gradients. Indeed, such a change 
could not be appreciated by our observations as those central H/3 
profiles are in fact unresolved. On the other hand, if we were to 
advocate the role of AGNs also where the H/3 surface-brightness 
profile is resolved and is slightly shallower than an r^^ law (e.g., 
NGC 2768, NGC 2974, NGC 4486) then projection effects would 
have to be invoked in addition to a radial decrement of the gas den- 
sity. Considering vertically extended structures may also help to 
explain some of the sharpest peaks for the value of the gas veloc- 
ity dispersions that are observed in our sample (e.g., NGC 3414, 
NGC 4278, NGC 5198). 

As a final item of this section, in Fig. [5] we draw the 
[N l]/H/3 vs. [O IIl]/H/3 diagnostic diagram for the ionised-gas 
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Figure 5. SAURON [Nl]/H/3 vs. [O III]/H/3 diagnostic diagram for galaxies with 1.4 Ghz radio cores. The different symbols con'espond to each of the 
galaxies with detected [N I] emission, and the size of the symbols indicate the radial range over which the average values of the [N I]/H/3 and [O in]/H/3 
ratios have been computed. The colored grids show the same MAPPINGS-111 models as in Fig[T] except for the AGN photoionisation models shown in the 
left panel, which do not include dust, con'espond to a single value of the power-law index of the AGN continuum of a = —1.2 and are shown also for solar 
values of the gas metallicity. The range of the dimensionless ionising parameter of the dust-free AGN models is also somewhat different, including values of 
log U = -3.3, -3.0, -2.6, -2.3, -2.0, -1.6. 



emission observed in the central 3" of the SAURON galaxies 
with 1.4 Ghz radio-continuum cores and [Nl] emission. The 
emission in the central regions of these objects is generally 
consistent with both photoionisation from a central AGN with 
a very low ionisation parameter and with excitation from fast 
shocks. If we disregard the latter possibility for the moment (see 
^3.41 . we note that only models including gas o f super- solar 
metal l icity and a self-consistent tre a tment of dust jPopita et al.l 
I2002I ; |g roves. Dopita, & Sutherland! |2004|) provide a plausible 
description for our data. In particular, dust-free models with gas 
of solar metallicity (Fig. |5] left) can match just a fraction of our 
data (e.g., the very center of NGC 4486) and only with the hardest 
ionising radiation fields (i.e., oc i/^^'^). Super-solar values for 
the gas metallicity are required also in the case of dusty, radiation 
pressure-dominated photoionization (Fig. [5] right), since also these 
models would have to invoke the hardest ionising continuum in 
order to match our data in the case of gas of solar metallicity (not 
shown for the sake of clarity). The very nuclear regions of the 
objects in Fig. |5] (within 1", smallest symbols) appear consistent 
with dusty super-solar AGN models while adopting similarly low 
values of the ionisation parameter and a realistic range for the 
hardness of the AGN radiation field (from cx to u~^'^), 

in keeping with the behaviour of the SDSS LINER nuclei ( ij2.1| l. 
Figure |5] shows that significant radial variations in the average 
values of the [N l]/H/3 and [O IIl]/H/3 ratios are found only in 
NGC 2768, NGC 4278 and NGC 4486. As the average values of 
the line ratios do not run parallel to the dusty AGN models, such 
gradients are inconsistent with a simple radial change in the mean 
value of the ionising parameter (due to either a different mean 
distance or gas density) or of the gas density, suggesting instead 
a variation in the hardness of the AGN radiation field. This seems 
rather unlikely, however, in particular for NGC 4486 where the 
radial increase of the [N l]/H/9 ratio would imply a hardening of 



the AGN continuum hitting clouds at increasing distance from the 
centre. Therefore, excluding the possibility of special gas and dust 
geometries, we consider it more likely that other processes than 
AGN photoionisation are responsible for the gas excitation beyond 
the very centre of these obj ects. In fact, f or NG C 4486 the results 
of Sabraetal. (2003) and iDopita et all ( 1 19971) suggest a radial 
transition from nuclear AGN photoionisation to circumnuclear gas 
excited by fast shocks. 

To conclude, the presence of a 1.4 Ghz radio-continuum core 
suggests that a central AGN could be responsible for powering the 
ionised-gas emission in the very centre of 30% of our SAURON 
sample galaxies, at most. For AGN photoionisation to be consistent 
with the observed radial profiles of the [O IIl]/H/3 line ratio within 
the central 3", the central gas distribution has to have a radially 
decreasing density, and in the case of objects with resolved H/3 
emission profiles, also a vertically extended structure. Even so, the 
presence in few objects of considerable gradients in the value of 
[N l]/H/3 ratio suggests that other mechanisms than photoionisation 
by a central AGN must contribute to power the nebular emission 
outside the very centre. 



3.2 H ii-regions 

A few early-type galaxies in the SAURON sample show clear ev- 
idence of on-going star-formation. In particular, the values of the 
[O IIl]/H/3 ratio in NGC 3032, NGC 4526, NGC 4459 plunge to 
levels that are well below the lower-limit of the [O IIl]/H/3 ratio 
that is observed also in LINER-like emission regions and Seyfert 
nuclei, which is around [O IIl]/H/3 ~ 0.3 or log([0 IIl]/H/3)= 
-0.5 (see, e.g.. iKewlev et al.l l2006h . As noticed already in Pa- 
per V, these objects are also characterised by extremely regu- 



8 Sarzi et al. 



lar dust morphologies and gas kinematics, in particular that of 
the H/3 line, suggesting a relaxed gas distribution and dynam- 
ics. Similarly relaxed structures are observed also in NGC 524 
and NGC 5838, which too show low [O IIl]/H/3 values although 
never down to the low level that is observed in the previous 
galaxies. That star formation may be occurring in NGC 3032, 
NGC 4459, and NGC 4526 is supported also by the finding of 
massive disks of molecular g as ( Combes . Young. & Bureaull2007l : 
lYoung. Bureau & Cappellari |2007 ) and the detect ion of emis- 
sion f rom polycyclic aromatic hydrocarbons (PAHs IShapiro et al.l 
l2009t) . 

Yet, even in these relatively simple objects we observe strong 
gradients in the value of the [O IIl]/H/3 ratio, which underscore 
the presence of other sources of ionisation than OB-stars. In this 
respect it is useful to observe in Fig. [6] the behaviour in the 
[N l]/H/3 diagnostic diagram of the SAURON data for NGC 3032 
and NGC 4526, the only two star-forming early-type galaxies with 
detected [N I] emission. For comparison, in Fig.[6]we also show the 
same SAURON measureme nts for the Sa galaxy NGC 4314 (from 
iFalcon-Barroso et alj2006h , which hosts a central star-forming ring 
and where the nebular emission is stronger, in terms of equivalent 
width of the lines, than in NGC 3032 and NGC 4526. Figure |6] 
confirms dramatically the presence of star formation in NGC 3032 
and NGC 4526 as it shows that the SAURON measurements in 
these two objects display the same characteristic diagonal trend 
in the [N l]/H/9 vs. [O IIl]/H/9 diagram that is observed also in 
SDSS galaxies with star-bursting nuclei (Fig.[T]l and in the ring of 
NGC 4314. Furthermore, Fig.|6]confirms the impression conveyed 
by the maps for the equivalent width of the H/3 emission (Fig. 4b 
of Paper V) that, as in the case of NGC 4314, star formation in 
NGC 3032 and NGC 4526 occurs primarly in a ring. Indeed, in 
Fig. |6] the SAURON data reveal a striking similarity between the 
line-ratio gradients observed in NGC 3032, NGC 4526 and the Sa 
NGC 4314, whereby in all three galaxies the emission-line proper- 
ties move closer to the predictions of the dusty AGN or shock mod- 
els as we consider regions away from the star-forming ring both 
towards the centre and the outer regions of the galaxy. In fact, in 
NGC3032 the regions outside the ring seem to approach the AGN 
and shock grids from a different direction than the points inside the 
ring, which would suggest either that toward the centre AGN pho- 
toionisation is more important than shock excitation, or that shocks 
may occur under different conditions at t he opposite edges of the 
ring. A similar trend was observed also bv lMazzuca et al.l ( l2006h in 
the Sa galaxy NGC 7742 using standard BPT diagrams. 

In Fig. |6]we also overplot on the SAURON maps the contour 
corresponding to the 1.4 Ghz radio-continuum emission measured 
by the FIRST survey. Consistent with synchrotron radiation from 
relativistic electrons and free-free emission from H II regions, the 
radio continuum is resolved and matches well the distribution of 
the nebular emission . Similarly to the case of normal disk galax- 
ies jKennicuttl ri983h . the 1.4 Ghz emission in NGC 3032 and 
NGC 4526 is dominated by the non-thermal synchrotron compo- 
nent, given that the thermal free-free emission that would corre- 
spond to our H/3 fluxes falls short to explain the ob served fl u x den- 
sity values. In fact, following th e prescription of ICondonl ( Il992h 
and ICapi an & Deharven 3 ( Il986l) to convert the H/3 emission into 
free-free emission, the observed 1.4 Ghz flux densities appear 
remarkably consistent with our SAURON emission-line measure- 
ments (accounting for beam smearing) if we adopt a 10% fraction 
for the thermal component t o the total 1.4 G hz fluxes, as found 
on average in disk galaxies jKennicuttl [T983h . That also in early- 
type galaxies the 1 .4 Ghz continuum traces mostly the synchrotron 



emission may explain why in Fig.[6]the radio contours do not ap- 
pear to follow in detail the H/3 emission observed in NGC 3032 
and in the Sa NGC 4314, in particular where the H/3 fluxes peak. 
Judging the agreement between radio and nebular emission is more 
difficult in NGC 4526, due to the higher inclination of its gaseous 
disk. It is instructive to note also the different impact of redden- 
ing by dust in NGC 3032 and NGC 4526 due to inclination effects, 
which explains why the latter galaxy shows brighter radio contin- 
uum flux densities, which are unaffected by dust, despite displaying 
fainter H/3 emission. 

No diffuse radio emission is detected in the rings of NGC 4459 
and NGC 524, due to the fact that the nebular emission in these 
objects is significantly fainter than that observed in NGC 3032 
and NGC 4526, and the corresponding radio emission has es- 
caped detection in the FIRST images. Given the typical detec- 
tion threshold of 0.45 mjy/beam in the FIRST survey, we ex- 
pect to detect radio emission associated to star-formation in the 
FIRST images only if the average surface brightness of the H/3 
line within 99.73% of the FIRST beam (FWHM=5'.'4) exceeds 
~ 0.9 X 10~^''erg cm~^s~^arcsec~^, including also the contri- 
bution of synchrotron emission. This level of surface brightness for 
the H/3 line (conveniently shown in Fig. 4a of Paper V in light- 
green colours) is never achieved in the ring regions of NGC 524 
and NGC 4459. On the other hand, the H/3 and [O III] flux values 
in NGC 5838 are comparable to what is measured in NGC 4526, 
suggesting that the unresolved 1.4 Ghz emission shown in Fig. |3] 
could arise from the compact gaseous disk of NGC 5838 as well 
as from a central AGN, in agreement also with the fact that this 
object sits on the far-IR to rad io-continuum emission relation 
dCombes. Young. & Bureaull2007l) . 

Could photoionisation by OB-stars be responsible also for the 
nebular emission that is observed in galaxies with less relaxed dust 
and gas structures or log([0 IIl]/H/3) values greater than —0.5? Star 
formation is unlikely to be the source of excitation for clouds with 
log([0 IIl]/H/3) > 0.5, corresponding to starburst models with low 
metallicity and large q values (Fig.[Tll. Star formation episodes char- 
acterised by subsolar metallicity values are indeed extremely rare 
(because gas is quickly enriched during such events), whereas the 
ionisation parameter of extragalactic H Il-regions is typically found 
to range only between 10^ < g < 10** (Dopitaetal. 2000). For 
more intermediate values of the [O IIl]/H/3 ratio, the possibility that 
OB-stars could be powering the observed emission is particularly 
compelling when there is evidence for recent star formation, for in- 
stance from the analys is of optica l spectra (Paper VI and VIII) or of 
near-UV colours (e.g.. |jeong et aL 20 02), or in the presence of mas - 
sive reservoirs of molecular gas jCombes. Young. & Bureaul2007h . 
In this respect, cases like NGC 3489 and NGC 3156 are particularly 
intriguing (but see also Wi.5\ . Both of the previous conditions are 
met in these objects, where the gas kinematics is also fairly relaxed 
with small values of the gas velocity dispersion o-gas. Yet, despite 
these similarities with the star-forming objects, the [O IIl]/H/3 ra- 
tio in NGC 3489 and NGC 3156 is consistently > 1, peaking to 
remarkably high values (~10) in places well outside the nucleus, 
which rules out H II regions as the main source of nebular emission. 

For these and other objects with evidence of recent star 
formation or reservoirs of molecular gas we can use the VLA 
FIRST radio-continuum data to further test whether the observed 
nebular fluxes are entirely associated to star formation activity. 
For instance, if the nebular emission observed in NGC 3489 was 
powered by OB-stars, following the previous prescription to derive 
1.4 Ghz fluxes from the observed H/3 emission, we should see a 
well resolved radio structure, as in the case of the star-forming 
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Figure 6. Two examples of lenticular galaxies with extremely regular dusty disks and gas kinematics where star formation is on-going, shown together to the 
Sa galaxy NGC 4314 which hosts a well-known circumnuclear star-forming ring. The top and middle panels show maps for the H/3 and [N I]AA5198, 5200 
emission. Superimposed to these maps are 1 .4 Ghz radio continuum contours (red lines) from the FIRST- VLA images illustrating the resolved character of the 
synchrotron and free-free emission. As in the Paper V and for the SAURON maps in following figures, the white contours follow the stellar surface brightness, 
the scale is in arcsec and the side box show the range of the plotted quantity, in this case the flux of the H/3 and [N I] lines, in logarithmic scale. The [N I] 
doublet is detected only where [O III] and H/3 emission is the strongest, which allows to trace these regions in [N I]/H/3 diagnostic diagrams shown in the 
lower panels. Similar to the case of the most intense starburst observed in SDSS data (Fig.[T), all three galaxies share a distinct diagonal trend in the [N I]/H/3 
diagnostic, with the central and outermost regions (shown in darker and lighter points, respectively) departing the most from the starburst model grids and 
falling closer to the AGN or shock grids. 



galaxy NGC 4526 that indeed displays similar values for the 
H/3 fluxes. This is not observed, since NGC 3489 presents a 
rather weak and only marginally resolved core ( ii3.lt Similarly, in 
NGC 3156 the central H/3 emission is sufficiently intense that if it 
were associated to star formation it would lead to a detectable radio 
core, but this is not observed. The non detection of 1.4 Ghz radio- 
continuum features excludes the possibility that nebular emission 
is powered solely by OB-stars also in the nuclei of NGC 2685, 
NGC 4150 and NGC 7332, which have young stellar populations, 
relatively bright central H/3 emission and, except for NGC 7332, 
also molecular gas. Among the objects with old stellar ages, as 
estimated from the SAURON absorption-line measurements, it is 
worth considering the cas es of NGC 2768, which has a r obust 
single-dish CO detection jCombes. Young. & Bureaull2007l) , and 
of NGC 2974, where the near-UV colours sug gest the presence o f 
young stars in the centre and in an outer ring jjeong et"al]|2007h . 
Although both galaxies display strong radio-continuum cores 



that are undoubtably dominated by AGN activity, the observed 
H/3 fluxes are sufficiently strong and extended (exceeding the 
~ 0.9 X 10~^''ergcm~^s~^arcsec~^ threshold. Fig. \^ that if 
they were entirely powered by OB-star the corresponing free-free 
and synchrothron emission from star-forming regions would be 
capable of inducing noticeable and extended features in the radio- 
continuum maps. That in NGC 2768 the radio core is unresolved 
therefore suggests that star-formation cannot be principal the 
source of ionization in the central regions. Conversely, the fact 
that the radial profile for the 1.4Ghz continuum in NGC 2974 
is resolved does not necessarly imply that the H/3 emission is 
powered by OB-stars. In fact, this possibility will be ruled out in 
^3.41 by means of the SAURON [Nl]/H/3 diagnostic, suggesting 
that the extended nature of the radio core of NGC 2974 may be 
due to the presence of an unresolved jet. Nonwithstanding the lack 
of evidence for nuclear star formation, we note that NGC 2768 
and NGC 2974 are among the brightest IRAS far-infrared sources 



10 Sarzi et al. 



in our sample (together with NGC 524, NGC 3032, NGC 4459, 
NGC 4526 and NGC 5838) suggesting that star-formation occurs 
primarily outside the centre in these two galaxies. 

To conclude, extremely low values of the [O IIl]/H/3 ratio 
leave no doubt that OB-stars are powering most of the nebular 
emission observed in the 6% of the SAURON sample galaxies 
(3/48), with regular dust morphologies and relaxed gas kinematics 
strongly suggesting that star formation could be responsible for the 
gas ionisation in up to 10% of the sample (5/48). The presence of 
young stellar populations, molecular gas, and the detection of PAHs 
features and strong FIR fluxes further support the case for on-going 
star formation in these objects. For two of them, NGC 3032 and 
NGC 4526, the SAURON [N l]/H/3- diagnostic diagram reveal a dra- 
matic similarity with the star formation activity in the central rings 
of the early-type disk galaxies, showing also the transition to re- 
gions outside the ring where other excitation mechanisms appear to 
dominate, such as AGN photoionisation towards the centre. The de- 
tection in the FIRST data for NGC 3032, NGC 4526 and NGC 5838 
of radio-continuum emission consistent with synchrotron radiation 
from relativistic electrons and free-free emission from H Il-region 
further supports the case for on-going star-formation in these ob- 
jects, whereas the non-detection of extended 1.4 Ghz emission in 
NGC 524 and NGC 4459 can be ascribed to a lower star-formation 
activity and the limited sensitivity of the FIRST survey. The FIRST 
data also suggest that in other objects with young stellar popula- 
tions and molecular reservoirs only a fraction of the H/3 fluxes can 
be attributed to H 1 1 -regions. 



3.3 Diffuse Evolved Stellar Sources 

Apart from very young stars, several kinds of evolved stars or stel- 
lar remnants emit significant amounts of ionising photons, such 
as post-asymptotic giant branch stars (pAGB), extreme horizontal- 
branch stars and X-ray binaries. Although none of these objects 
produce an ionising continuum as powerful as that of a single O- 
star, their importance as photoionisation sources lies in their large 
numbers and potential ubiquity, in particular since early-type galax- 
ies are made mostly of old stellar populations. In fact, the sig- 
nificance of pAGB star s as ionising candidates was first demon- 
strated bv iBinette et al. 

I fl994), who showed that when taken to- 
gether such objects produce a sufficient Lyman continuum to ac- 
count for the Hg luminosity of early-type galaxies. Shortly after- 
wards, |Macchetto|eran^ finding of a correlation 
(alreadv hinted at bv lPhillips et"ar]|l986l) between the Ha-l-[N II] 
luminosity of early-type galaxies and their integrated B-band lumi- 
nosity within the typical extent of the nebular emission, thus adding 
to the hypothesis that the principal sources of ionising photons are 
to be found within the bulk of the stellar population. 



3.3.1 The Connection Between Nebular and Stellar Emission 

Following iMacchetto et al.l we computed for each of our sam- 
ple galaxies the average distance from the center of the detected 
emission-line regions, and integrated the flux of both the stellar 
continuum and of the H/3 recombination line observed over all the 
SAURON spectra within such a radius. After rescaling for galactic 
distance, the left panel of Fig. [7] shows a nearly linear correlation 
between the luminosity of the H/3 line and that of the stellar contin- 
uum within the average extent of the nebular em ission, consistent 
with the narrow-band results of IMacchetto et al] 



Yet, the true strength of the link between the nebular and stel- 
lar emission in early-type galaxies is dramatically revealed when 
the total H/3 luminosity is compared with the luminosity of the stel- 
lar continuum if this is integrated only where the nebular emission 
is detected. As the central panel of Fig.|7]illustrates, in this case the 
correlation between the values of the H/3 luminosity Lh/3 and con- 
tinuum luminosity Lcont is much tighteiQ. Such a relation is nearly, 
although not quite, linear (with a power-law slope a — 1.30±0.08), 
meaning that the equivalent width of the integrated H/3 emission is 
almost constant across our sample. The tightness of Lcont~Lu0 re- 
lation (with a 0.15 dex scatter) is quite remarkable considering that 
in addition to a diffuse stellar sub-population other sources of ion- 
isation may contribute to the total flux of the H/3 recombination 
line. In fact, the only obvious outlier in Fig. [7] is NGC 3032, the 
galaxy displaying by far the most intense star-forming activity in 
our sample. The possible contribution of a central AGN does not 
affect much the Lcont-Lup relation, as its power-law slope and 
scatter do not change significantly (a = 1.20 ± 0.05, 0.13 dex 
scatter) when the central 3" are excluded ( i]3.1t . The tightness of 
this correlation is not a distance artifact either, since it persists also 
when the flux of the H/3 emission and of the stellar continuum are 
compared (with a — 1.31 ± 0.08), as shown in the right panel of 
Fig. [7] A direct comparison of the nebular and continuum fluxes 
also dispels the possibility that the observed trend is due to a sensi- 
tivity bias. As discussed in Paper V, our detection thresholds trans- 
late into a rather narrow range for the minimum equivalent width 
of the emission lines that we can measure, from 0.2A to just 0.07 A 
for H/3, so that in practice also the flux of barely detectable H/3 
lines scales with the strength of the underlying stellar continuum. 
Reassuringly, the right panel of Fig. |7] shows that the integrated 
fluxes corresponding to such detection limits not only lie comfort- 
ably below the total measured values of the H/3 flux but also that 
they define a shallower correlation with the continuum fluxes than 
found in the case of the observed H/3 fluxes. 

The tight correlation between the integrated luminosity of the 
nebular emission and that of the stellar continuum stems in re- 
ality from the fact that within most of the galaxies in our sam- 
ple the radial profile of the ionised-gas emission follows very 
closely the stellar surface brightness distribution. Figure[8]provides 
three specific examples (in the case of NGC 2974, NGC 3414 and 
NGC 4150) that vividly illustrates how well the nebular fluxes can 
follow the stellar continuum in early-type galaxies, or alternatively, 
of how the equivalent width of H/3 can take remarkably similar val- 
ues at very different distances from the galactic centre correspond- 
ing to very different continuum levels of the stellar surface bright- 
ness (in contrast to the case of NGC 3032, also shown in Fig.^. In 
fact, across our entire sample the equivalent width of the H/3 line 
fluctuates on average by just 32% around the mean value it takes 
in any given galaxy, excluding objects with obvious star formation 
(NGC 3032, NGC 4459, NGC 4526; Sg^J and while disregarding 
the central 3" of all galaxies to avoid any possible AGN contamina- 
tion ( i]3.1t . This is conveyed also by Fig.|9]where the H/3 and con- 
tinuum fluxes contributed by each bin of our sample galaxies (with 
the same exceptions) are plotted against each other, after rescaling 
the H/3 and continuum fluxes observed in each galaxy to common 
average values. A power-law fit to such a stacked distribution indi- 



^ With the benefit of hindsig ht a tigh t correlation can also be observed in 
the measurements of Macche tto et al.l if only the objects with diffuse and 
extended emission (dubbed "DE") are considered. 



The SAURON project -XVI 11 




Figure 7. The connection between the integrated luminosity of the nebular emission and that of the stellar continuum in early-type galaxies. Left: Correlation 
between the luminosity of the stellar continuum Lcont and that of the H/3 emission Lh/3 observed in the SAURON spectra when both quantities are integrated 
within the average radius of the emission-hne regions detected in each of our sample galaxies. The solid and dashed line show the best fitting power-law 
fit to the observed distiibution and associated Icr scatter, respectively. Middle: Same as left but now the total H/3 luminosity is juxtaposed to the stellar 
luminosity integrated only where the nebular H/3 emission is detected. A much tighter correlation emerges in this case. Right: Same as rniddle but now 
showing the correlation between the integrated nebular and stellar fluxes. A tight correlation persists, indicating that the Lcont-iH/3 relation is not driven by 
the distance distribution of our sample galaxies. The light grey points show the integrated H/3 detection limits in each galaxy. Their position and alignment 
differs significantly from that of our integrated H/3 fluxes, dispelling also any doubts that the observed correlations are due to the limited range of our equivalent 
width sensitivity. In all panels the position of NGC 3032, our most significant outlier, is indicated, and can be explained in terms of intense star-formation 
activity. 




Figure 8. Examples of the local connection between the flux of the nebular emission and stellar suiface brightness in early-type galaxies. Maps for the 
equivalent width of the H/3 recombination line {top panels) display remarkably similar values for this quantity across most of the SAURON field of view, 
reflecting how closely the flux of the H/3 line follows that in the stellar continuum (lower panels). Such a correlation breaks down in the central regions of our 
sample galaxies (shown with darker points in the lower panels) where AGN photoionisation can be important ( il3.1K or in objects where OB-stars are the main 
source of ionisation (e.g. in NGC 3032, rightmost panels) 



cates a nearly linear relation (a = 0.88 ±0.01) with 40% scatter 
(0.14 dex). 

3.3.2 Explaining the Nebular to Stellar Connection 

Does a nearly constant equivalent width for the recombination lines 
necessarily imply a diffuse stellar photoionising source? After all, 
the nebular flux and the stellar surface brightness, the two elements 
entering the equivalent width definition, relate to two different inte- 



grals of the stellar luminosity density. Whereas the surface bright- 
ness is a line-of-sight integration of the luminosity density, the re- 
combination flux is connected to the total Lyman continuum emit- 
ted by the putative ionising stellar sub-population, which reaches 
the gas from all directions within the galaxy. 

The simplest model for the radial profile of the equivalent 
width of H/3 can be computed in the case of a spherical galactic 
geometry and assuming that a constant fraction of the ionising radi- 
ation from the evolved stellar sources is absorbed and subsequently 
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Figure 9. Correlation between the flux of the steUar continuum and the 
flux of the H/3 hne (left panel) observed in each of the bin of our sample 
galaxies where emission is detected, except for the central 3" and excluding 
the objects with obvious star formation. For each of the plotted galaxies, the 
continuum and nebular flux values have been rescaled to the same average 
value, of 2.9 X 10~" and 6.4 X 10" ■^^ergs ^cm ^, respectively. The 
right panel shows the same correlation, but for the detection limits on the 
H/3 line, normalised by the same factor used for the observed H/3 values. 



re-emitted across the galaxy by the ionised-gas, where the latter 
is also taken to reside in a planar configuration. Under the previ- 
ous assumptions the equivalent width of the recombination lines at 
a given projected distance R from the center will be directly pro- 
portional to the ratio of total stellar radiation -F(ii) to the stellar 
surface brightness at the same point, since also the stellar 

ionising continuum reaching the gas clouds should scale with the 
stellar radiation F(R) (for a set stellar population). For any given 
intrinsic luminosity density profile L(r) the integral that delivers 
FiJC) is, in spherical coordinates, 



dr 
J a 



^^r^ sin Odd 



(1) 



where x — \/ R'^ + — 2rR cos 9. Substituting cos 9 = jj,, dfi = 
— sin 9 d9 the previous triple integral can be reduced to 
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Adopting a 7-model for the intrinsic luminosity density of early- 
type galaxies (Deh neiilll993h 



L(r) = 



(3-7)M 



47r 



(4) 



rT(r + a)4-T 

where M is the total luminosity (or stellar mass) of the galaxy and 
a is a scaling radius (which relates to R^ through 7), the integral 
given by Eq.[3]becomes 
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This can be numerically evaluated by breaking it in two pieces: 
from r = Q ^ R substituting r = Ru, and from r = R ~* 
00 substituting r = R/u. The 7-models also allow to s olve for 
the s urface brightness profile (see Appendix B of lOehnenl 

1 1993 ). which can be varied to include fairly shallow models (7 = 1, 
corresponding to Hernquist 1990 profile) as well as rather steep 
profiles (7 = 2, corresponding to a Jaffe 1983 profile). 

With both F{R) and E(i?) at hand, Fig.[Tn]shows that accord- 
ing to our simple model the equivalent width of the recombination 
lines should increase towards the outskirts of early-type galaxies. 




Figure 10. Predicted radial profile for ratio of the ionising radiation F(R) 
from evolved stellar sources diffuse throughout the galaxy to the surface 
stellar brightness I{R). Assuming that a constant fraction of F{R) is re- 
processed as H/3 emission, the plotted F/I profile provides a prediction 
for the radial variation for the equivalent width of H/3. The model further 
assumes a spherical geometry and three different ki nd of 7-models for the 



intrinsic luminosity density L (r), approximatin g a ide Vaucouleun 



Jaffe 



1948, 



1983 



7 = 1.5) or con'esponding to Hernquist fl99d, 7 = 1), and 
7 = 2) surface-brightness profiles. Independent of the chosen L{r), the 
equivalent width of the recombination lines should increase towards the 
outskirts of early-type galaxies according to this simple model. 



with little dependence on the choice of the intrinsic L{r) profile. 
Between ~10% and ~70% of one effective radius R^, which on 
average correspond to 3" and the radius containing 90% of the 
emission in our sample galaxies, the model predicts a ~80% rise 
in the equivalent width of H/3. Although a 0.25 dex increase may 
be accomodated by the average trend observed across our sample 
(Fig.[9l(, more often than not the equivalent width of the H/3 line 
tend to remain constant or even decrease with radius, in particu- 
lar if we focus on the objects with the most extended and diffuse 
H/3 emission (e.g.. Fig. [8]( and exclude those where the impact of 
other sources of ionisation than evolved stars may be particularly 
important (e.g., the filaments of NGC 4486). 

Can such a discrepancy be remedied by relaxing the simple 
premises of our model? An important aspect to keep in mind is that 
not all the ionising photons emitted by the evolved stellar sources 
may reach the gas due to interstellar absorption, so that only the 
stellar radiation emitted within the mean-free path I of the photons 
should be considered. This will lead to a steeper F{R) profile and a 
shallower F{R) /S(_R) gradient, in better agreement with our data. 
Yet, as ; ^ and F{R) L{r) ■ I the F/I ratio will eventually 
decrease as fast R~^, again in disagreement with our observations. 
As regards our assumption that a constant fraction of the ionising 
photons is absorbed across the galaxy, this depends on how closely 
the radial profiles of the gas column density and of ionising ra- 
diation follow each other, and on whether the filling factor of the 
gas clouds varies. The column density of a steady disk increases 
towar ds the center of a real i stic ga lactic potential well, as the mod- 
els of lKawata. Cen.&Hol j2007l) illustrate. The work of Kawata 
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Figure 11. lonisation balance for post-AGB stars, low-mass X-ray binaries and extreme horizontal branch stars, within the region containing 90% of the 
nebular emission of our sample galaxies. Darker shades of grey indicate galaxies with an higher area coverage of the ionised-gas, where the ionising flux 
from the stellar sources should be re-proccessed more efficiently. Objects where O- and B-stars have been previously estiblished ( ^3.2t as the main source of 
ionisation (NGC 3032, NGC 4459, NGC 4526) are excluded from these diagrams. 



et al. (2007) also shows that outside the very central regions the 
Toomre parameter Q is rather constant (see their Fig. 2), so that 
the gas surface density profile Igas{R) should scale almost exactly 
with the epicyclic frequency k (given that Q oc n/Igaa)- Using the 
7-models to compute n one learns that the gas surface density de- 
creases always more gently than F{R) with radius, and even more 
so when compared to L{r). For a fixed filling factor it is therefore 
likely that an increasing number of ionising photons escape the 
galaxy as we approach the centre, whereas towards the outskirts 
of the gaseous disk an increasing fraction of the gas may remain 
neutral. The combination of a small mean-free path I and of a flat- 
tening of the H/3 flux profile towards the centre (as it departs from 
the F{R) ~ L{r) profile) could thus provide a relatively simple 
explanation for the general behaviour observed across our sample 
galaxies of a nearly constant equivalent width of H/3. Alternatively, 
if the ionising photons travel far enough in the galaxy that F{R) 
is not significantly affected by I, another plausible avenue to ex- 
plain our data is to consider the possibility that the ionising sources 
are more concentrated than the bulk of the stellar populations. This 
may be the case if the ionising stars are also responsible for the 
far-UV flux (~ 1500A) observed in early-type galaxies. The far- 
UV light displays indeed a steeper surface-brightness profile than 
observed at optical wavelengths, following closely als o the metal- 
licity gradient traced by the Mgb line-strength index jjeong et"ail 
l2009t) . 



3.3.3 lonisation Balance 

In the absence of a definite model for the tight connection between 
the values of the stellar and nebular luminosity of early-type galax- 
ies, we turn our attention to the energy budget. Specifically, we ask 
whether the various ionising stellar candidates can provide a suffi- 
cient number of ionising photons to explain the observed H/3 fluxes. 
In doing so we return to the use of integrated flux measurements, 
since within large apertures the stellar luminosity estimated from 
the integrated flux approximates well the total luminosity within 
the galaxy, and therefore provides a basis for estimating the total 

ionising stellar radiation. 

As mentioned above, iMacchetto et alj ( 1 19961) already carried 
out the ionisation-budget exercise in the case of pAGB stars. Using 



the predictions o f lBinetteetanh994h for the specific ionising pho- 
ton luminosity from such stars (that is, per unit mass of the entire 
stellar population) and assuming t hat all the ion ising radiation is 
re-processed into nebular emission, IMacchetto e t al. could convert 
the stellar luminosity within the average radius of the emission- 
line regions into a limiting value Lna.pAGB for the Hq luminosity 
that could be po ssibly powered by pAGB stars. Here we revisit 
'Macchetto et al. experiment, using the stellar population models 
of Yi et al. ( 1999) to compute the Lyman continuum from pAGB 
stars and 2MASS K-band images for our sample galaxies to trace 
the stellar mass within the radius containing 90% of the nebular 
emission. Consistent with pr evious w ork, the specific number of 
pAGB stars in the models of lYi et al] is pretty much constant be- 
tween 4 and 11 Gyr, which brackets the typical age of the stel- 
lar population dominating the stellar mass of our objects. Assum- 
ing case B recombination the specific Lyman luminosity of pAGB 
stars in these models would deliver 0.007 ergs"^ in Ha photons 
for each K-band ergs^^, if all ionising photons are intercepted 
and re-processed. Under these circumntances, FigurefTTIshows that 
pAGB stars still provide an abundant number of ionising photons 
to explain the ionised-gas emission observed in our sample galax- 
ies. That the predicted limiting values I/Ha,pAGB far exceed the 
observed values Lua is encouraging considering that for many ob- 
jects the emission-line regions cover only a fraction of the area used 
to compute such upper-limits (shown with a whiter shade of grey), 
so that almost certainly a good fraction of the pAGB ionising radia- 
tion escapes the galaxy. On the other hand, we need to keep in mind 
that the present stellar population models may be grossly overes- 
timating the specific number of pAGB stars, and therefore their 
ability to power the ionised-gas emission in early type galaxies. 
The recent near- and far-UV obs ervations of M32 with the Hubble 
Space Telescope by iBrown et al.i (2008) suggests indeed a dearth 
of pAGB stars compared to standard predictions for their numbers. 

Turning our attention to other stellar candidates, we now con- 
sider the role low-mass X-ray binaries (LMXBs). The exquisite 
spatial resolution of the Chandra X-ray observatory has recently 
allowed to detect LMXBs in nearby ellipticals, enabling the deriva- 
tion of their luminosity function and, most important in the present 
context, their total X-r ay luminosity (in the 0.3-8 keV range; e.g. 
iKim & Fabbiandl2004h . According to lKim & Fabbianol the LMXB 
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population within the optical radius of their host galaxies can de- 
liver from a few times 10'^^ to several lO^^ergs"^. The X-ray 
spectral energy distribution of LMXBs in early- type galaxies has 
been found to be well ch aracterised by a simple power-l aw with 
photon index F — 1.56 jlrwin. Athev. & Bregman|[2003h . corre- 
sponding to an energy index a = 0.56 (for Fi, oc v~"). For a 
power-law spectrum with an energy index a — 0.1 to 0.9 the Ha 
luminosity produced by reprocessing the X-ray emission can range 
from 1/9 to 1/2 of the X-ray luminosity (Ho 2008). Thus, to a 
first approximation, the overall LMXBs population of early-type 
galaxies could provide enough ionising photons to power the ob- 
served nebular emission, given that the Hq luminosity of the latter 
ranges between a few times 10'^* and several lO'^^ergs"^ (see Pa- 
per V). Encourag ed by this first check we used the relation found by 
iKim & Fabbiand between the K-band luminosity Lk of early-type 
galaxies and the X-ray luminosity Lx,lmxbs of their LMXB pop- 
ulation to more accurately estimate for each of our sample galaxies 
the limiting Ha luminosity Lhq.lmxbs that could be powered by 
these stellar remnants. Upon such a closer inspection, and assuming 
that 18% of the LMXBs X-ray luminosity converts into Ha emis- 
sion, Fig.ll llshows that LMXBs do not live up to our expectations 
as ionising sources in early-type galaxies. In particular, Fig.ll llre- 
veals that although the Lhc«,lmxbs values are on average consistent 
with the observed Ha luminosities, in agreement with our previous 
estimates, in practice none of our sample galaxies could display 
emission powered solely by LMXBs if the gas coverage is taken 
into account. In fact, in all objects with I/Hq.lmxbs ^ Lua in 
Fig. [TT| most the LMXBs ionising radiation is likely to escape the 
galaxy because of a patchy gas distribution, whereas when the gas 
coverage is high the LMXB population seems unable to power the 
nebular emission. 

Finally we explore the case of extreme horizontal-branch stars 
(EHBs). Such low-mass core Helium burning stars only appear in 
significant numbers in old stellar systems, contributing in partic- 
ular to the far-UV flux in populations with ages greater than 9 
Gyr. For this reason, EHBs have been suggested as the dominant 
sourc e for the UV-upturn obs erved in many early-type galaxies 
(e.g., lGreggio & Renzin llggOl) and as indicators of the oldest stel- 
lar populations jYi et alJ ll999Y Although EHBs are much weaker 
as ionising sources than pAGB stars, their contribution to the Ly- 
man continuum becomes nonetheless significant in stellar popula- 
tions older than 10 Gyr, eventually overtaking pAGBs as the main 
source of ionising photons when the population turns 12 Gyr old. 
In order to trace the contribution of EHBs in our sample galaxies 
we used GALEX measurements for the far-UV flux (1341-1809 A), 
extracted within the area containing 90% of the nebular emission as 
in the case of our 2MASS measurements. Contrary to the case of 
the near-UV passband (centred at ~2310A), the far-UV measure- 
ments are fairly insensitive to the contribution of younger stellar 
components (i.e. from few 100 Myr to 1 Gyr), whereas pAGB stars 
do not contribute to either of these passbands. To transform the 
GALEX far-UV fluxes into a limi ting Hq l uminosity from EHBs, 
i'Hof.EHBs, we used the models of lYi et aU to integrate the Lyman 
continuum due only to stars either on horizontal-branch or that have 
only recently evolved from it, which includes also AGB-manque 
stars. We adopted models for a 1 1-Gyr-old stellar population, which 
reproduce the UV-to-optical colour observed in prototypical UV- 
upturn galaxies such as NGC 4552. Metallicity does not impact on 
our results as this parameter only produces differences in the mod- 
els in UV regions (~ ll OOA) out side the far-UV passband. In this 
framework the models of iYi et al] yield 0.0295 erg in Ha pho- 
tons for each far-UV ergs~^. As in the case of LMXBs, Fig. II II 



shows that EHBs per se are unlikely to generally power the nebular 
emission observed in early-type galaxies. 



3.3.4 Summary on the Role of Evolved Stellar Sources 

To conclude, we have added to early suggestions that the main 
source of ionisation for the gas of early-type galaxies is to be found 
among the bulk of their stellar population. We have found a very 
tight correlation between the total values of the nebular and stellar 
luminosity in our sample galaxies, when considering only the re- 
gions where the ionised-gas emission is detected. Furthermore, we 
have shown that such a global trend stems actually from the fact 
that the ionised-gas emission of early-type galaxies follows very 
closely their stellar surface brightness. Although we have not been 
able to conclusively explain how such a constant profile for the 
equivalent width of the nebular emission can be produced by stellar 
photoionisation, through simple modelling and reasoning we have 
shown the direction for future investigations without invoking em- 
barassing fine tuning of key parameters such as the mean free path 
of the ionising photons. Overall, pAGBs remain the favourite ion- 
ising candidates to explain the gas emission of early-type galaxies, 
which points to a need for further investigations into the present 
discrepancy between observations and models as regards their total 
numbers. The role of LMXBs would seem marginal in this context, 
unless future observations reveal a population of super-soft X-ray 
sources, which would ionise more efficiently the gas surrounding 
them and induce a very specific nebular spectrum ( Rappaport et al.l 
1 19941) . In principle, the relative importance of pAGBs, LMXBs and 
EHBs could be further investigated through a careful emission- 
line diagnostic, since this analysis can trace the specific shape 
of the ionising continuum of the different sources. For instance, 
EHBs have a much softer Lyman continuum compared to pAGBs, 
which should make it more difficult for EHBs to produce extended 
partially-ionised regions. Unfortunately a diagnostic analysis based 
on the [N I] emission is not possible for the diffuse component dis- 
cussed here, since the [N I] lines are detected almost exclusively in 
the central regions of our galaxies, in the presence of star formation, 
or in the filamentary structures of the most massive of our sample 
galaxies where the interaction with the hot ISM is the most likely 
source of ionisation (i ]3.6l l. Notwithstanding the absence of diffuse 
[N I] emission, the behaviour of the [O IIl]/H/3 ratio is still useful in 
this context. Specifically, we note that strong gradients and peculiar 
features are not expected to arise easily in the presence of a uniform 
radiation field from the same kind of source across the galaxy, un- 
less to consider contrived structures of higher dust concentration or 
enhanced gas number density that could allow for local fluctuation 
of the ionising parameter q. The finding of remarkable [O IIl]/H/3 
structures suggests therefore the role of additional ionising mecha- 
nisms, such as those discussed in the next sections. 



3.4 Shocks 

3. 4. 1 Fast Shocks 

Shocks offer a natural explanation for some of the most prominent 
and coherent structures that are observed across our sample in the 
maps for the equivalent width of the H/3 and [O III] lines and for 
the [O IIl]/H/3 ratio, in particular for morphologies like spiral arms, 
oval and integral-sign shaped structures (e.g., Fig. I12t. Such fea- 
tures can indeed be induced by the presence of stellar bars or weak 
triaxial perturbation of the stellar potential, which can funnel the 
gas into preferential streams where fast (above 100 km s^^) shocks 
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between gas clouds may occur. Yet, the presence of equivalent- 
width structures per se does not guarantee that shocks are pow- 
ering the gas emission - other diffuse sources (such as pAGB stars) 
may be still ionising the gas as this accumulates in spiral arms and 
similar features, where the observed emission would be stronger 
only due to a higher column density. On the other hand, the case 
for shock excitation becomes much stronger when the features in 
the [O III] or H/3 emission are accompanied by similar structures 
in the [O IIl]/H/3 ratio. Similarly, we should expect a correlation 
between the strength of the recombination H/3 emission and the 
velo city dispers ion of the lines in the presence of fast shocks (see, 
e.g. jDopita & S utherland 1995), because the ionising flux powered 
by the latter scales with the shock velocity and because unresolved 
velocity gradients could in principle be traced by an enhanced ve- 
locity dispersion. 

Associated [O IIl]/H/3 structures are more often observed 
for integral-sign emission features (where the [O in]/H/3 distribu- 
tion can be oddly asymmetric, e.g., NGC 4262, NGC 4278 and 
NGC 4546), and only in NGC 3414 in the case of spiral arms. 
As regards the gas velocity dispersion CTgas, this quantity does 
not seem to generally trace notable equivalent-width features, ex- 
cept for NGC 4546 where o-gas appears to display larger values 
along the locus of the main integral-sign structure. That remarkable 
equivalent-width structures are not always accompanied by simi- 
larly obvious [O IIl]/H/3 and (jgas features (as Fig. [12] illustrates 
for (Tgas) suggests that shocks may indeed not always be the main 
source of ionisation in these regions. On the other hand, the shock 
geometry may not lead to significant broadening (as set by the spec- 
tral resolution of the SAURON data) of the emission lines along the 
line of sight, for instance if the shocks occur in regions that are 
much smaller than the physical scales that we can resolve (typically 
■~100pc) and if they are surrounded by a more relaxed ionised-gas 
medium. In this respect, it is interesting that all galaxies featuring 
notable spiral arms or integral-sign structures display also extended 
central velocity dispersion peaks (see also Fig. I12t . Such central 
(Tgas gradients could mark the point where the previous limitations 
are overcome and the kinematic signature of the shocks can be read, 
if the shock fronts wind up and get closer together as the gas is fun- 
neled towards the center. This certainly appears to be the case for 
NGC 2974, where Emsellem, Goudfrooij, & Ferruit (2003) showed 
with high-resolution integral-field spectroscopic observations that 
the spiral-arm morphology observed in the SAURON gas maps ex- 
tends well towards the centre. Emsellem et al. could also resolve 
better the central Ugas gradient in the SAURON data, finding (jgas 
peaks on the inner sides of the central spiral arms that further sug- 
gest the presence of shocks, although yet not necessarily that these 
can power the ionised-gas emission. 

When [N I] emission is detected in objects with notable spiral 
or integral-sign features, the corresponding [N l]-diagnostic analy- 
sis is always consistent with shock-excitation. The [N l]-diagnostic 
is seldom conclusive, however, since most of the [N I] emission 
originates close to the centre of our sample galaxies (see e.g., the 
cases of NGC 2974 and NGC 4278 in Fig. [Hi where the [N I]- 
diagnostic diagrams indicate that also AGN activity could be pow- 
ering the gas emission. In fact, if shocks are responsible for the gas 
excitation in these regions, the grids of Fig.[T2]suggest exceedingly 
high shock velocities between 300 and 500 kms~^, which can 
hardly be explained throug h gravitational motions. For instance, the 
gas-dynamical models of lEmsellem. Goudfrooij. & FerruitI j2003l) 
for the inner regions of NGC 2974 suggest gas streamlines crossing 
at an angle of just 45° (see their Fig. 21), which even assuming gas 
motions at nearly circular speed (up to ~300 km within the 
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Figure 12. Two examples of elliptical galaxies where the presence of 
spiral-arms, ovals, or integral-sign morphologies for the ionised-gas distri- 
bution highlights the impact of non-axisymmetric perturbation of the grav- 
itational potential on the gas motions, possibly leading to shocks between 
gas clouds. Emission from the [N I] doublet (second row of panels) is de- 
tected where the [0 III] (upper panels) and H/3 emission is the strongest, 
allowing for a limited emission-line diagnostic analysis (lower panels). The 
colored grids show the same MAPPINGS-III models as in Fig[T] whereas 
the grey grids display also the predictions for shock models with a gaseous 
precursor ahead of the shock front. These last grids are for models with the 
same values of the shock velocity Vs and the magnetic parameter b as for 
the shock models without precursor H II regions (see the caption of Fig. [T). 
The [N I]/H/3 diagnostic indicates that where the [N I] lines ai'e detected in 
NGC 2974 and NGC 4278 the observed emission is consistent with both 
AGN photoinisation and shock excitation, although in the latter case the 
required values for Vb would range between 300 and 500 km s~^, which 
can hai'dly be attained thi'ough gravitational motions (see text). Note that 
increasing the electron density of the schock models (n^ = 1) would only 
shift the grid to the left as the [N I] lines collisionally de-exites more easily, 
which would imply even higher Vs values to power the observed emission 
with shocks. Maps for the gas velocity dispersion Cgas (third row of pan- 
els) are also shown to illustrate how equivalent width structures (top) are not 
always accompanied by obvious crgas features, as expected if o-gas would 
trace Vs . 
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central 2-3" according to the models of Paper IV) would not allow 
to reach the required shock velocities. It thus seems unlikely that 
shocks are powering the nebular emission observed in the central 
[N l]-emitting regions of these two galaxies. 

Despite the previous shortcomings shocks could still have a 
role in powering the gas emission observed in integral-sign fea- 
tures when such structures are accompanied by notable asymmetric 
[O IIl]/H/3 distributions, as is the case for NGC 4262, NGC 4278 
and NGC 4546 (see Fig. 4 of Paper V). Excluding AGN photoion- 
isation for such extended structures and star-formation, which usu- 
ally proceed in relaxed morphologies such as rings or disks, ex- 
plaining such disparate [O IIl]/H/3 values through photo-ionisation 
by a diffuse source such as pAGB stars would require different val- 
ues for the ionising parameter q at the opposite sides of a galaxy and 
therefore having either different values for the gas number density 
or for the concentration of dust. Maintaining such different physical 
conditions over time would seem unlikely at first, which leaves us 
to consider the role of shocks. One way to explain an asymmetric 
distribution for the [O IIl]/H/3 ratio is to suppose that shock ioni- 
sation proceeds under different conditions ahead and behind of the 
shock fronts and that dust precludes us from receiving most of the 
nebular emission from either side of such fronts, which will nat- 
urally occur when looking at one en d or the other of the galaxy . 
In fac t, shock models such as those o f iDopita & SutherlandlTim 
1 19961) feature the possibility of adding a gaseous precursor ahead of 
the shocks, which would naturally occur in gas-rich galaxies. Such 
precursor H II regions are characterised by high ionisation, which 
makes them efficient [O III] emitters and therefore good candidates 
as the principal source of nebular emission along the line of sight 
where the highest [O IIl]/H/3 values are observed. Fi sure [T2l shows 
the predictions also of the shock models with a gaseous precursor, 
which indeed can explain the observed values for the [O IIl]/H/3 ra- 
tio that typically range between ~ 2 and 4. Yet, also in the presence 
of a precursor the model grids imply rather large high Va values be- 
tween 200 and 300km s^^, which admittedly could still be difficult 
to achieve, in particular in objects with well organised gas motions 
like NGC 4278. 

In fact, our premise to discard the role of a diffuse source of 
ionisation in order to explain such asymmetric [O IIl]/H/3 struc- 
tures, that an uneven dust distribution would be unlikely, is actually 
unfound ed according t o Spitzer observations. Based on IRAC im- 
ages, in Ishapiro et alj we uncover features in maps for the 8pim 
non-stellar emission excess that neatly correspond to the regions 
in NGC 4278 where the [O IIl]/H/3 ratio falls at its lowest value. 
Such a 8/im flux excess is usually associated across our sample 
to star-forming regions, and comes from ionised-PAH emission. 
In NGC 4278 there is no evidence for on-going star formation, 
however, in particular given the absence of a substantial molecu- 
lar reservoio The 8/im excess observed here may be due to pre- 
dominantly neutral PAHs that emit more significantly in the 11.3 
and 12.7//m bands, similarly to other early-type galaxies with PAH 
e mission and no CO detection (se e, e.g., the case of NGC 2974 
in iKaneda. Onaka. & Sakonll2005h . Such an asymmetric concen- 
tration of dust would induce a local hardening of the Lyman con- 
tinuum from evolved stellar sources in the galaxy and would re- 
duce the rate at which the ionising photons hit the gas clouds, thus 



^ The single-dish detection of ICombes. Young. & Bureau! )2007l) in 
NGC 4278 is only marginal compai'ed to other galaxies showing evidence 
for on-going (e.g., NGC 3032, NGC 4459 and NGC 4526) or recent (e.g., 
NGC 3156 and NGC3489, see 333) star formation. 



leading to a smaller ionising paramet er q (for a discussion on dust 
effects see lShields & Kennicutll 19951 and references therein). Both 
effects would produce lower values for the [O IIl]/H/3 ratio in re- 
gions with the highest dust covering. 

Finally, shocks would also spring to mind as the most likely 
source of ionisation for the emission observed in galaxies like 
NGC 3156 and NGC 3489, when glancing at the extreme range 
of values for the [OlIl]/H/3 ratio (from ~1 up to 10) that is 
found across the entire SAURON field of view in these objects 
(see Fig. |13b . Yet, the high [O IIl]/H/3 values observed in these 
objects are inconsistent with shock excitation if we consider that 
NGC 3156 and NGC 3489 are among the least massive galaxies 
in our sample. The circular velocities in these objects peak at just 
~150 km s~^, which would allow only for a restricted range of 
shock velocities (say Va ^ 200 km s^^) and thus only for a rather 
limited range of values for the [O IIl]/H/3 ratio (up to 2, see Fig.ll2t. 
irrespective of the presence of a gaseous precursor. Furhermore, 
NGC 3156 and NGC 3489 display fairly regular gas kinematics that 
indicates the presence of a relaxed gaseous system, which rules out 
the possibility of strong non-circular motions to circumnvent the 
previous impasse with shock excitation. 



3.4.2 Slow Shocks 

The dynamical deadlocks encount ered in the previous su bsection 
may suggest that slow shocks (e.g.. lShull & McKedl 19791) . instead 
of fast shocks, could be important. Slow shocks could be driven by 
ram pressure into gas clouds as they move through the hot halo gas, 
and have the appealing property of producing relatively strong low- 
ionisation lines such as [O I] or [N I], and values of the [N l]/H/3 
ratio close to what we observe. Even assuming that all early-type 
galaxies can retain a halo of hot. X-ray emitting gas, whic h is far 
from clear to be the case (e.g.. lMathews & Brighentill2003l but see 
also i]3.6t , a scenario involving slow shocks would run into two 
sorts of complications. 

First, the values of the [O IIl]/H/3 ratio predicted by the 
IShuU & McKe3 models are quite sensitive upon the shock veloc- 
ity, so that in practice slow shocks would be consistent with our 
data only for a rather restricted range of shock velocities, between 
85 and 95 km s~^. Such a condition seems unlikely to be gener- 
ally satisfied, however, considering the wide range of velocities at 
which the gas clouds may be orbiting within our sample galaxies. 

Second, if slow shocks were generally occurring in early-type 
galaxies, they could lead to a profile for the equivalent width of the 
recombination lines that would be inconsistent with the flat profiles 
that we observe, and even more than what we simply estimated in 
the case of photoionisation by old stars ( i]3.3.2| l. If the local density 
of the hot gas sca les like the square root of the stellar density (as 
in giant ellipticals. lTrinchieri. Fabbiano. & CanizareslT986) we can 
compute the radial trend that the emission from slow shock would 
have, using the same simple assumptions and models of i]3.3.2| and 
considering that the line emission from such shocks scales like the 
cube of the orbital velocity of the clouds (presumably close to cir- 
cular) and the density of the hot halo gas. Within the same region 
between 10% to 90% of an effective radius that we considered in 
i]3.3.2| the equivalent width of emission arising from slow shocks 
would increase markedly a factor 2.5, which is not supported by ob- 
servations. Furthermore, the equivalent width of lines arising from 
slow shocks would drop to zero towards the centre, posing another 
inconsistency with the data in the absence of additional central 
source such as an AGN. 
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Figure 13. Two prototypical examples of early-type galaxies with intense 
nebular emission characterised by large fluctuations for the [O III]/H/3 ra- 
tio, which also peaks to relatively high values throughout extended regions. 
From top to bottom, the panels show maps for the equivalent width of the 
[O III] line, for the [O III]/H/3 ratio and for the gas velocity dispersion, re- 
spectively. Objects displaying such highly-ionised extended regions (HI- 
ERs) are not very massive and show invariably evidence for a recent star- 
formation episode. Both the intensity of the emission and the large val- 
ues for the [O III]/H/3 ratio observed in HIERs can be attained consider- 
ing as ionising sources the pAGB stars associated to the recently formed 
stellar population (see ^3.5K whereas the large fluctuations in the observed 
[O ni]/H/3 values is explained considering that star formation may still be 
occurring in regions with very low [O III]/H/3. The relaxed character of 
the gas kinematics, as traced by the extremely small values of the gas ve- 
locity dispersion, suggests that the conditions for star formation are still 
favourable in the specific case of these two objects. 

3.4.3 Summary on the Role of Shocks 

To conclude, shocks, either fast or slow, can hardly be regarded as 
an important source of ionisation for the diffuse nebular emission 
observed in early-type galaxies. In particular, we have discussed 
how fast shocks could contribute to power the ionised-gas emis- 
sion in regions where the gas is funneled into spiral and integral- 
sign structures by the non-axisymmetric perturbations of the gravi- 
tational potential. Yet, even there the role of diffuse stellar sources 
may be as, if not more, important than shocks, in particular when 
the spiral and integral-sign features in maps for the [O III] or H/3 
emission are not accompanied by similar [O IIl]/H/3 structures. In 
the central regions of our sample galaxies shocks could contribute 



to explain the high values for the gas velocity dispersion, but lim- 
ited [N I] diagnostic analysis combined with simple dynamical ar- 
guments appear to dismiss shocks as the main source of gas ex- 
citation. Arguments for the importance of shocks, which included 
also the presence of a gaseous precursor and of dust, were put for- 
ward also to explain the peculiar asymmetric [O IIl]/H/9 distribu- 
tion that is often associated to integral-sign emission structures, but 
ultimately run into a similar dynamical dead lock. In these cases 
Spitzer observations reveal an uneven distribution of dust, which 
can help reconcile the observed [O IIl]/H/3 character of the nebu- 
lar emission with simple photoionisation by diffuse stellar sources. 
Finally, we have also considered the case for shock excitation in 
low-mass galaxies with high-ionisation extended regions, and even- 
tually ruled out shocks as the main source of ionisation. In the next 
section we will consider an alternative explanation for such extreme 
cases. 



3.5 Post Starbursts and High-Ionisation Extended Regions 

An intriguing subset of the early-type galaxies surveyed in the 
course of the SAURON survey consists of relatively small objects 
(with (Je ^ 120km s~^) displaying high-ionisation extended re- 
gions (HIERs), where the [O IIl]/H/3 ratio reaches some of the 
most extreme values observed in our sample with great varia- 
tions across the SAURON field of view (Fig. [13}. NGC 3156 and 
NGC 3489, which we already considered in i]3.2| and i]3.4| repre- 
sent only the most remarkable examples of such a kind of early- 
type galaxies - HIERs are also found in NGC 7332 and NGC 7457 
(see Fig 4 of Paper V). A common trait that these low-mass ob- 
jects share is the presence throughout the galaxy of fairly strong 
H/3 absorption features (Paper VI), which is indicative of a re- 
cent star-formation episode. In fact, star formation in NGC 3156 
and NGC 3489 may have not yet come to a complete halt, in 
particular in those relatively small regions where very low val- 
ues of the [O IIl]/H/3 ratio are observed. Both PAH and CO reser- 
voirs haTC_beeii^etectedJn_^^ (see Shapiro et al. 2009 
and ICombes. Young. & Bureaull2007l . respectively), although the 
amount of molecular material found in NGC 3156 and NGC 3489 
(Mh2 = 0.2,0.1 X 10*Mq) is significantly less than what found 
in the presently star-forming galaxies NGC 3032, NGC 4459 and 
NGC 4526 (with AfH2 = 2.5, 1.7, 3.7 x 10** M©). Additionally, the 
gas kinematics in NGC 3156 and NGC 3489 trace a remarkably 
cold dynamical system (the values of the gas velocity dispersion 
(jgas observed in NGC 3156 and NGC 3489 among the smallest 
measured in our sample, ~50 km s~^) which could still favour 
the formation of stars. All these facts suggest that these objects 
are post-starbursting sy stems, perhaps linked to "E-l-A" galaxiefl 
( lDressler&Gunnlll983h . 

Is the high-ionisation character of the nebular emission from 
these objects linked to t heir recent star-formation history? In the 
context of AGN activity, iTaniguchi. Shiova. & Muravamal ( |200(]|) 
considered the role of post-starbursting stellar systems to explain 
the LINER and Seyfert emission of nearby galactic nuclei. Namely, 
Taniguchi et al. recognised that as the most massive and recently 
formed stars start leaving the asymptotic giant branch (~ 10**yrs 
after the starburst), the corresponding planetary nebulae nuclei be- 
come an important and sufficiently hot source of ionising photons. 



^ The very nuclear regions, within the few tens of pc, of NGC 3489 already 
display the classical features of an E+A spectrum ISarzi et alj2005l) . 
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Whereas this scenario is fairly limited in explaining AG Ns , in par- 
ticula r since most galactic n uclei are invariably old (e.g. JSarzi et alj 
I2OO5L but see also lHoll2008 ). photo-ionisation from the pAGB pop- 
ulation associated to a recent starburst is a very likely ionising 
mechanism for the HIER emission observed in low-mass early-type 
galaxies. 

According to Taniguchi et al., to power with pAGB stars 
from massive progenitors (between 3 to 6 Mq) the Ha luminos- 
ity of NGC 3156 and NGC 3489, which is around ~ 10^''ergs~\ 
these objects should have experienced a starburst involving approx- 
imately ~ IO^Mq of gas, this is consistent with our estimates for 
the amount of stars that should have recently formed in these galax- 
ies (corresponding to 10% of the total mass) in order for the values 
of their stellar and dynamical mass-to-light ratio to agree with each 
other as is the case for more quiescent early-type galaxies in the 
SAURON sample (Paper IV). 

Another attractive feature of such a post-starburst photo- 
ionisation scenario is that it could naturally account for the 
extreme values of the [O IIl]/H/3 ratio observed in early-type 
galaxies with HIERs. The pAGB population associated with the 
recent starburst is presumably still confined to the galactic plane 
where star formation took place and where the remaining gas is 
still orbiting, thus bringing the gas much closer to its photoionising 
sources than is the case for the pAGB stars related to the older 
bulge population ( i]3.3t . In turn this situation would translate into 
an higher ionisation parameter q and therefore to larger values for 
the [O IIl]/H/3 ratio. We note also that due to their proximity to the 
gas, the Lyman continuum of the pAGB stars associated to a recent 
starburst could be reprocessed as nebular emission more efficiently 
than for other stellar sources throughout the bulge. This would help 
explaining also the remarkable strength of the nebular emission in 
objects like NGC 3156 and NGC 3489, where the equivalent width 
for the integrated H/3 and [O III] lines reaches values of 0.5A and 
1.5 A, respectively. 

To conclude, in approximately 10% of the early-type galaxies 
in the SAURON sample the nebular emission is most likely pow- 
ered by pAGB stars associated to a recently formed population 
of stars. The intense and highly-ionised character of the observed 
emission in these objects can be explained by the vicinity of the 
gas to its ionising stellar sources, which are presumably still con- 
fined to the galactic plane where they formed. The nature of the 
nebular activity observed in the post-starbursting systems of the 
SAURON sample could also have bearings for our understanding 
of the "E-h A" phenomenon. "E -l-A" galaxies are extremely rare ob- 
jects (e.g., iKavirai et al.1l2007h that display strong Balmer absorp- 
tion lines but no nebular emission, suggestive of a recent but al- 
ready exhausted, perhaps even quenched, star-formation episode. 
Yet, "E-l-A" are not totally dev oid of gas since they often still dis- 
play abundant H I reservoirs ( BuvleetalJl2006) . which, according 
to our picture, could be re-ionised just a few ~ lO^yrs after the 
starburst. This would leave only a very short time after the death 
of OB-stars for the galaxy to appear quiescent while developing 
A-star type features, which could contribute to explain why so few 
post-starbursting galaxies are observed in their "E-i-A" phase. 



3.6 Interaction with the Hot Interstellar Medium 

Historically, it was perhaps the finding of optical emission in and 
around the most massive galaxies of X-ray selected clusters that 
first drew considerabl e attention to the ionised - gas component 
of early-type galaxies tau. Cowie. & Wandfigsl : iHeckman et al.l 
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Figure 14. X-ray properties of SAURON early-type galaxies. The total 
blue-band luminosity Lb is compared to the X-ray luminosity Lx, nor- 
malised to Lb- Symbols are colour-coded according to the obse rved de- 
gree of rotational support as traced by the Ajj pai'ameter of Emsell em et alj 
i2007), with green or blue symbols showing slowly- or non-rotating ob- 
jects with Afl < 0.1. Lx values (circles) or upper limits (small cir- 
cles with downward arr ows) are ROSAT and Einstein m easurements 
from the compilation of lO' Sullivan, Forbes. & PonmanI j2 00l|), excep t for 
NGC524 and NGC58 13. which were taken from IPellegrina j200a) and 
iBohringer et alj fcOOOl) . respectively. These data cover three-quarters of the 
SAURON early-type sample, with the remaining objects being typically low- 
luminosity galaxies with an average log Lb ~ 9.9. The horizontal dashed 
line shows the normalised X-ray lum inosity expected from stellar sources 
<0' Sullivan. Forbes. & Ponm^ll200l! ). Only luminous galaxies can retain 
massive X-ray halos, in particular if they are slowly- or non-rotating ob- 
jects. 



Il989l) . The discovery of such a warm medium was indeed regarded 
as a major success of the cooling flow theory (Fabian 1994), al- 
though it soon became apparent that the amount of ionised gas 
found in the optical filaments of X -ray clusters exceeded by far the 
rates predicted by cooling flows jHeckman et al.|[T989h . Further- 
more, the presence of dust in such filaments posed an additional 
problem to the cooling flow scenario since gas cooling down from 
a KeV plasma should not be dusty and since dust particles ought 
to be quickly destroyed b y X-ray spu ttering in environments per- 
meated by hot ions (e.g.. lDonahue & Voit 1993). The connection 
between the X-ray and optical emission in and around the most 
massive early-type galaxies remained an established fact, however, 
which led to the suggestion that the nebular emission arises in fact 
from material that has been recently accreted and that is being 
excited by t he hot medium itself, e ither thro ugh thermal electron 
conduction jSparks. Macchetto. & Golombekl Il989l : Ide Jong et al.l 
1 19901) or by the radiati on that the hot gas emits as it cools down 
i Voit&Donahuelll990h . Sparks et al. (1989) showed that the en- 
ergy flux provided from heat condu ction is su fficient to explain 
the nebular emission, whereas Donah ue & VoitI (1991) found that 
the X-ray gas radiation could explain also the typical line ratios 
obse rved in the optical filaments of cooling flows clusters (see 
also iFerland et al.ll2009l for a recent investigation). Yet, the most 
spectacular confirmation of the connection between the hot and 
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Figure 15. Warm versus hot gas emission in the four brightest and slowly- or non-rotating galaxies in the SAURON galaxies, which all show also prominent 
dust lanes. The top panels show SAURON maps for the equivalent width of the H/3 line, in A and on a logarithmic scale, whereas the lower panels show 
Voronoi-binned maps for the X-ray emission observed with Chandra, in the 0.3-5.0 keV band. Except for only a few regions, the nebular emission is generally 
associ ated to X-ray emitting features. For NGC 4486 and NGC 5846, this connection was already discussed bv .Sparks et al.. 1 .2004.) and ,Trinchieri & Goudfrooii 
<2002l) . respectively. 



warm phases of the ISM came when X-ray images of high spatial 
resolution revealed a striking spatial coincidence between regions 
displaying optical and X-ray emission, first in the c ase of M87 



jYoung, Wilson, & MundellH2002l: ISparks et al.ll2004h and then in 
NGC 5846 jTrinchieri & GoudfrooiilboO^) ! 

In light of these findings we turn our attention to the inter- 
action between the warm and hot phases of the ISM, as the last 
ionising mechanism that we will consider in order to explain the 
nebular emission observed across our own sample of early-type 
galaxies. Unsurprisingly, it is only the most luminous galaxies in 
the SAURON sample that can retain a massive halo of hot, X-ray 
emitting gas, although the degree of rotational support also appears 
to play a part in determining the presence of an X-ray atmosphere. 
Figure [14] vividly illustrates this point by showing how only the 
brightest galaxies featu re values of the X-ray luminosity (m ostly 
from the compilation of lO' Sullivan. Forbes. & Ponmanll200ll) that 
well exceed what is expected from an unresolved population of 
LMXBs, in particular if considering slowly- or non-rotating galax- 
ies (green or blue circles). All but two (NGC 4552, NGC 5982) 
of such massive systems display prominent dust features (see the 
HST maps of Paper V), and luckily Chandra images are available 
for these dusty giant galaxies (NGC 4374, NGC 4486, NGC 5813, 
NGC 5846) in order for us to investigate the connection between 
the X-ray and optical emission. 

Figure [Tsl compares the Chandra X-ray images for the four 



brightest slowly- or non-rotating galaxies in the SAURON sam- 
ple with maps for their H/3 emission from Paper V. The Chan- 
dra images show integrated fluxes in the 0.3 - 5.0 keV energy 
band and have been spatially r ebinned through Voronoi tessela- 
tions jCappellari & CoDinll2003h to reach a minimum of 10 counts 
per bin. Except in only a few regions, the ionised-gas emission is 
usually associated to X-ray emitting features, whereas the converse 
does not always hold. Thus, with NGC 4374 and NGC 5813 we 
bring two additional examples besides NGC 4486 and NGC 5846 
of objects displaying such a rema rkable coincidence b etween the 
hot and warm phases of the ISM. 'S parks et alj ( |2004|) could also 
show that the X-ray filaments of NGC 4486 have a lower temper- 
ature than the surrounding medium, supporting the idea that heat 
is being transferred from the hot to the warm gas. Unfortunately 
the Chandra images for the other objects in Fig. [15] are not suffi- 
ciently deep to allow the same kind of analysis, since binning up 
to reach the minimum signal quality for a sensible spectral anal- 
ysis (usually 50 counts) would leave only a few patches. Yet, our 
integral-field data reveal another characteristic of the ionised-gas 
emission in such dusty giant slowly- or non-rotating galaxies that 
sets them aside from the rest of the SAURON sample and which re- 
inforces the case for gas excitation through the interaction with the 
hot medium. 

Figure [T6] shows indeed how the optical filaments of such 
dusty giant ellipticals display rather intermediate values for 
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Figure 16. Median values and Icr confidence limits for the [O III]/H/3 ratio 
observed in the SAURON galaxies with clearly detected and extended H/3 
and [O III] emission, plotted in alphabetical order against the galaxy name. 
The brightest, dusty and slowly- or non-rotating ellipticals are shown with 
green symbols, objects with dusty disks and on-going star formation are 
shown in blue, and the remainder of this subsample is denoted in black. 
The horizontal dotted and dashed lines show the average and Icr bound- 
ary for the [O III]/H/3 ratios observed in the galaxies shown by the black 
points. The tendency for the dusty giant ellipticals to show lower values for 
[O III]/H/3 ratio than for the rest of the quiescent galaxies in the SAURON 
sample suggests, in light of the presence in these objects of a massive halo 
of X-ray emitting gas, that the interaction between the hot and warm phases 
of the ISM has an impact on the ionisation of the nebular plasma. 



the [Oin]/H/3 ratio (~1-1.5), which lie below the average 
[O IIl]/H/3 values observed in the other quiescent early-type 
galaxies in the SAURON sample. This low-ionisation pattern is 
reminiscent of the low [O IIl]/H/3 values that are observed in the 
optical filaments in and around the brightest cluster galaxies (e.g. 
ISabra. Shields. & Filippenko 2000: Hatch et al.l2o6^ . Yet, in such 
enviroiraients the [O IIl]/H/3 ratio falls to much lower values than 
we observe in our sample g alaxies when optic al and X-ray emis- 
sion coincide. For instance, iHatch et al] ( |2006|) hardly detect any 
[O III] emission in the filaments that extend well beyond the optical 
boundaries of NGC1275, the brightest member of the Perseus 
cluster, whereas in the inner region of the same ionised-gas nebula 
ISabra. Shields. & Filippenkol ( l200(]h measure [O III]/H/? ~ 0.3-1, 
still below the values that we find in the four galaxie s of Fig 1151 
Rather than considering varying Oxygen abundances, iHatch et al] 
interpret this behaviour as an ionisation gradient driven by a 
progressive intensification of the excitation mechanism toward the 
optical regions of NGC1275. We agree with this conclusion, and 
further suggest that the intermediate [O IIl]/H/3 values observed in 
the optical regions of our sample galaxies where hot and warm gas 
appear to interact are in fact caused by the contribution of stellar 
photoionisation, in addition to the interaction with the hot ISM. 
There is indeed no reason to exclude in dusty giant ellipticals the 
presence of a Lyman continuum from old stars, which per se would 



lead to the larger [O IIl]/H/3 values that are observed in the rest of 
the quiescent early-type galaxies of the SAURON sample. In fact, 
we note that for the optical filaments associated to X-ray features 
the radial profile of the H/3 flux tends to follow rather closely the 
stellar surface brightness, further suggesting the importance of 
stellar photoionisation. 

To conclude, in the brightest and slowly- or non-rotating early- 
type galaxies the presence of a massive halo of hot. X-ray emitting 
gas imprints a distinct low-ionisation character to the emission-line 
regions beyond the direct influence of a central AGN. This line- 
ratio pattern is achieved thanks to the energy input provided by the 
hot gas as it interacts, presumably through heat conduction, with 
the warm gas, which is still subject to photoionisation from old 
stellar sources like in other early-type galaxies. As reported in pre- 
vious work, also in the dusty giant galaxies of the SAURON sample 
a remarkable spatial coincidence between optical and X-ray emit- 
ting features testify to the interaction between the hot and warm 
phases of the ISM. 



4 IMPLICATIONS FOR LARGE-SCALE SURVEYS 

If the high sensitivity of the SAURON survey has revealed the extent 
to which diffuse ionised-gas emission is found in nearby early-type 
galaxies (up to ~ 75%, Paper V), shallower but much larger spec- 
troscopic campaigns such as the Sloan Digital Sky Survey (SDSS) 
can help reveal through shear numbers important connections be- 
tween the nebular and stellar properties of early-type galaxies (e.g., 
Kauffmann et al. 2003; Graves et al. 2007; Schawinski et al. 200"^; 
Kewlev et alj2006l ; lKauffmann & Heckmanl2009il The nebular ac- 
tivity detected in the SDSS spectra of early-type galaxies is, ac- 
cording to standard BPT-diagnostics (see also iQ, generally either 
LINER-like or characterised by composite H Il/Seyfert-like emis- 
sion (often referred to as transition objects, or TOs), especially 
when early-type galaxies are selected o n the basis of their re d 
colour rather than by their morphology jSchawinski et alj|2007h . 
Such a kind of n ebular emission is most often taken as due to a 
central AGN (e.g., Kauffmann e t alj2003l ; lKauffmaim & HeckmanI 
^009; Kewley et al. 2006), although unlike for H II and Seyfert nu- 
clei there are a number of plausible ionising mechanisms that can 
power composite and LINER-like emission. This is particularly rel- 
evant if one considers that the fixed 3"-wide aperture of the SDSS 
spectra usually encompasses large Kpc-scale regions where diffuse 
emission such as that observed across our sample could contribute 
to the SDSS nebular fluxes. In fact, other sources than a central 
AGN may significantly contribute to power the central LINER-like 
emission of even much closer galaxies, and therefore over much 
smaller physical regions than in the case of SDSS spectra, typic ally 
on 10 pc scales (see, e.g., the case of the Palomar survey, IHoI 
l2008l) . 

We can use our previous results on the relative importance of 
the various ionising mechanisms at work in and across different 
kinds of early-type galaxies in order to better understand the nature 
of the nebular emission observed in early-type galaxies in large- 
scale surveys, addressing in particular the real extent of AGN activ- 
ity. Indeed, we can exploit the integral-field nature of the SAURON 
data to simulate how our sample galaxies would appear if observed 
through wide circular apertures, as is the case for distant early-type 
galaxies targeted by the SDSS. A good starting point for such a 
comp arison between SAURON and SDSS data is the MoSES sam- 
ple of lSchawinski et ^ ( l2007h . which includes morphologically se- 
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lected early-type galaxies with redshift values between z — 0.05 
and 2 — 0.1. Using the SDSS pipeline parameter devRad.r as a 
proxy for the effective radius R^, in 68% of the MoSES objects the 
3"-wide aperture of the SDSS spectra encompasses between 32% 
and 79% of 7?e, an area which is almost always covered by the 
SAURON observations of our local sample. At these galactic scales, 
nebular emission is detected in the MoSES sample in 18.5% of 
the cases, with values for the equivalent width of the [O III] emis- 
sion typically (68%) between 0.8A and 4.3A, extending most often 
to 2.9A and 2.4A if we consider only objects with composite or 
LlNER-like emission, respectively. Seyfert nuclei and star-forming 
objects indeed make up the majority of the MoSES galaxies with 
the most intense emission. 

Figure [T7] shows the radial profiles for the integrated flux and 
equivalent width of the [O III] emission over increasingly large 
circular apertures, for the 9 SAURON galaxies that display suffi- 
ciently intense emission to be detected at the typical sensitivity of 
the MoSES spectra. Out of these 9 objects one is a star-forming 
galaxy (NGC 3032) whereas the remaining 8 show integrated val- 
ues of the [O IIl]/H/3 ratio in the log([0 IIl]/H/3)= - 0.5 range, 
which is typical of LlNER-like emission (e.g.. Ho et al. 1997, TOs 
can extend to lower [O IIl]/H/3 values). These 8 SAURON galax- 
ies cover fairly well the range of equivalent width values observed 
for most of the MoSES objects with composite and LlNER-like 
emission (in particular the latter), but fail to extend in the most in- 
tense regime occupied more often by MoSES Seyfert nuclei and 
strongly starbursting galaxies. To some extent this is to be expected 
considering that Seyfert nuclei are exceedingly rare, making up 
only 1.5% of the MoSES sample, and given that Hll-nuclei are 
most often found in low-mass early-type galaxies, which are par- 
ticularly under-represented in the SAURON sample since this is a 
represe ntative but incomplete sample of the nearby early-type pop- 
ulation jde Zeeuw et ai] |2002l) . Yet, this may also suggest that the 
most intense regime of emission is where truly nuclear LINER 
and composite activity can be found in SDSS spectra. In fact, 
of the SAURON galaxies with detectable LlNER-like emission at 
the MoSES sensitivity level only half show radio or X-ray sign- 
posts of AGN activity ( ij3.1l ). Furthermore, even in these cases 
the normalised profiles for the integrated flux of the [O III] emis- 
sion (Fig. [T7] top panel) illustrate how the regions dominated by 
AGN activity (within 3", according to the conservative estimates 
of i]3.1t would contribute at most between 30% and 55% of the 
total [O III] emission detected within the MoSES apertures. Thus 
the SAURON galaxies that, if placed at the distance of the MoSES 
galaxies, would display the typical values for the equivalent width 
of [O III] and for the [O IIl]/H/3 ratio of the MoSES LlNER-like 
galaxies are in fact mostly, if not completely, dominated by diffuse 
emission powered by other sources other than a central AGN, most 
likely post-AGB stars ( ^3[ . 

At first glance, the fraction of SAURON objects (9 out of 48, 
or 19 ± 6%) with detectable integrated emission at the MoSES sen- 
sitivity level (hereafter referred simply as "detectable") would also 
seem consistent with the fraction of the MoSES sample displaying 
nebular emission (18.5%), even when the fraction of star-forming 
systems (1 out 48, or 2 ±2% compared to 4.3%) or of galaxies with 
composite or LlNER-like emission (8 out 48, or 17± 5% compared 
to 12.6%) are considered separately. Yet, it should be noted that the 
objects with composite emission in the MoSES sample come typ- 
ically (68%) with log([0 IIl]/H/3) values between -0.25 and 0.28, 
whereas the 8 "detectable" and non-starforming SAURON galaxies 
show log([0 IIl]/H/3) values between 0.29 and 0.43, a range that is 
much closer to that of the MoSES galaxies with LlNER-like emis- 
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Figure 17. Radial profiles for the integrated flux (normalised to the to- 
tal, top panel) and equivalent width of the [0 III] A5007 emission (lower 
panel) for the SAURON galaxies that display sufficiently intense emission 
to be detected at the typical s ensitivity of the SDSS sp ectra for the MoSES 
early-type galaxy sample of Schawinski et al.l j2007t) . The integrated flux 
and equivalent width values are obtained over circular apertures with radii 
expressed in units of the galaxy effective radius Re- In the lower panel, 
the rectangle drawn with solid lines indicates where 68% of the values for 
the physical size of the galactic regions that are encompassed by the 3"- 
wide aperture of the SDSS spectra are located, while also showing where 
68% of the values for the equivalent width of the [O III] emission detected 
in the MoSES spectra can be found. In the same rectangle, the dotted and 
dashed horizontal lines show the upper extent of the region containing 68% 
of the values for the equivalent witdh of the [O III] emission observed only 
in the MoSES galaxies displaying composite and LlNER-like emission, re- 
spectively. In the upper panel, the dot-dashed horizontal lines show, for the 
galaxies with Radio or X-ray signposts of nuclear activity, the most conser- 
vative contribution of the AGN to the total [O III] flux (that is, within 3"; 
see i]3.1t . whereas the vertical lines show again the typical size of the galac- 
tic regions observed by the MoSES spectra, which can be used to estimate 
the AGN contribution to the nebulai' emission measured in such SDSS data. 



sion (where log([0 IIl]/H/3) spans between 0.18 and 0.55). To put 
it in a more direct way, the fraction of the MoSES sample that dis- 
plays log([0 IIl]/H/3) values in the - 0.5 range, which is generally 
associated with LlNER-like emission, is just 8.3%, or only about 
half the fraction of the "detectable" SAURON galaxies with inte- 
grated log([0 IIl]/H/3) values in the same range. It is possible that 
the lack of "detectable" SAURON galaxies with integrated values 
for [O IIl]/H/3 ratio that are more typical of composite activity is 
due to small-number statistics and the incomplete character of the 
SAURON sample. Indeed, as for the objects displaying star forma- 
tion and Seyfert activity, also the MoSES galaxies with composite 
emission live preferentially in low mass early-type galaxies, which 
are particularly under-represented in the SAURON sample. Alterna- 
tively, it is possible that some of the MoSES galaxies would not 
be classified as early-type if observed at a much closer distance. At 
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least for a local point of view, a more quantitative analysis will have 
to wait for a complete integral-field survey of the nearby early-type 
galaxy population, which will be possible once the data from the 
ATLAS^° campaigilflare at hand. 

Notwithstanding the limitations of the present analysis, the 
picture we have drawn for SDSS galaxies with [O IIl]/H/3 ratios 
typical of LINER-like emission (of which we have an over- 
abundance of "detect able" SAURON examples) would reinforce 
the recent results of IStasinska et id] {2OO8), who advocate the 
importance of photo-ionisation by post-AGB stars in "retired" 
galaxies to explain the LINER population in the SDSS survey. 
Yet, we note here that 2 out of the 8 "detectable" SAURON 
objects with LINER-like values of the [O IIl]/H/9 ratio are most 
likely powered by post-AGB associated to a recent star-formation 
episode (NGC 3156 and NGC 3489, ^33) rather than to the 
older bulge population. This kind of sources could be particularly 
common in the SDSS LINER population, considering that such 
SAURON post-starbursting systems appear to be more easily 
"detectable" in SDSS spectra (50%, or 2 out of 4) than is the case 
of the other SAURON galaxies with diffuse LINER-like emission 
but uniformly old stellar populations (25%, or 6 out 25). In fact, 
such a post-starburst L INER-like ac t ivity could contribute to 
explain the findings of I Graves et alj ( |2007|) . that red-sequence 
SDSS galaxies with strong LINER-like emission show younger 
stellar populations, in particular in less-massive systems, such as 
NGC 3156 and NGC 3489. 

To conclude, our SAURON integral-field data suggest that in 
very few, if any, of the SDSS galaxies which display only mod- 
est values for the equivalent width of the [O III] line (less than 
~2.4A) and LINER-like values for the [O IIl]/H/3 ratio, the neb- 
ular emission is truly powered by AGN activity. Only the most 
intense, and rarer, manifestations of LINER nuclear activity can 
be detected against both the stellar background encompassed by 
the SDSS aperture and the diffuse emission observed in early-type 
galaxies, unless to consider a very small minority of the closest 
SDSS objects. For SDSS samples that share these equivalent-width 
and line-ratio characteristics, the [O III]A5007 line is thus problem- 
atic as a diagnostic of low levels of black-hole growth, and LINERs 
may represent a mixed population that does not uniformly trace 
low-Eddington ratio accretion. 



5 CONCLUSIONS 

Building on our previous investigation of the incidence, morphol- 
ogy and kinematics of the ionised gas in the early-type galaxies of 
the SAURON sample, in this paper we set out to address the question 
of what is powering the nebular emission observed in them. 

To constrain the possible sources of gas excitation we resorted 
to a variety of ancillary data ranging from radio to X-ray wave- 
lengths, drew from complementary information on the gas kine- 
matics, stellar populations and galactic potential from our SAURON 
data, and introduced a SAURON-specific diagnostic diagram based 
on the [N l]/H/3 line ratio as a gauge for the hardness of the ion- 
ising continuum. Using SDSS data for galaxies with strong emis- 
sion we have shown that the [O IIl]/H/3 vs [N l]/H/3 diagram is a 
powerful tool to isolate star-forming regions and to clearly sepa- 
rate Seyfert and LINER nuclear activity. To best use the [N l]/H/9 

* See also' http : //purl . org/atlas3d| 



diagnostic we re-extracted more reliable [N I] measurements from 
the SAURON data using a new array of model and empirical tem- 
plates based on the MILES stellar library. Such empirical templates 
come from fits to high-quality SAURON spectra extracted over re- 
gions in most cases devoid of emission lines and with values for the 
absorption-line strengths that cover well the range observed across 
our sample. 

The considerable range of values for the [O IIl]/H/3 line ratio 
found both across the SAURON sample and within single galax- 
ies has prompted us to explore the relative importance of a central 
AGN, star formation, old UV-bright stars, shocks and of the interac- 
tion with the hot phase of the ISM in powering the nebular emission 
from early-type galaxies. The role of these ionisation sources can 
be summarised as follow: 

pAGB Stars. A tight correlation between the stellar surface 
brightness and the flux of the H/3 recombination line throughout 
the vast majority of our sample galaxies points to a diffuse and 
old stellar source as the main contributor of ionising photons in 
early-type galaxies. Based on simple ionisation-balan ce arguments 
we have shown that since the initial suggestion of iBinette et al.l 
( 1994) pAGB stars remain the favourite candidate for powering the 
ionised-gas emission of early-type galaxies, which urges further in- 
vestigations into the present discrepancy between observations and 
models as regards their total numbers. 

AGN. The presence of a radio or X-ray core suggests the ad- 
ditional role of AGN photoionisation in approximately 30% of the 
galaxies in our sample, although the radial profile of the H/3 recom- 
bination line shows that a central source can be responsible for the 
observed emission at best only within a few hundred parsecs from 
the centre, corresponding to the innermost 2" — 3". 

Shocks. Fast shocks are unlikely to be an important source of 
ionisation for the extended nebular emission observed in early-type 
galaxies, even where spiral or integral-sign gas structures are ac- 
companied by a similar [O IIl]/H/3 feature. Indeed, the shock ve- 
locities required to power the nebular emission can hardly be at- 
tained in the potential of our sample galaxies, in particular given 
the relatively relaxed gas kinematics that is observed in them. 

OB-stars. Extremely low values of the [O IIl]/H/9 ratio leave 
no doubt that OB-stars are powering most of the nebular emission 
observed in 6% of the SAURON sample galaxies (3/48), whereas 
the presence of regular dust morphologies, relaxed gas kinematics, 
young stellar populations, molecular gas, PAH features and strong 
FIR fluxes further support the case for on-going star formation in up 
to 10% of the sample (5/48). In these objects star formation appears 
confined mostly to circumnuclear regions, which in two cases is 
confirmed by [N l]/H/3 diagnostic diagrams that highlight the role 
of other sources of ionisation outside of the ring regions, for in- 
stance of a central AGN towards the centre. Yet, star formation may 
have proceeded also outside these circumnuclear regions, in partic- 
ular in the least massive galaxies like NGC 3032 where younger 
stellar populations are found throughout the entire SAURON field 
of view. 

Post-starburst pAGB Stars. For another 10% of the SAURON 
sample, consisting too of low-mass objects, the intense and highly- 
ionised emission relates also to a recent and spatially extended star- 
formation episode, although the principal ionising sources are no 
longer OB-stars but the pAGB population associated to the younger 
stellar subcomponent. 

Interaction with the Hot Phase. At the opposite end of the 
mass spectrum, we find that in the brightest and slowly- or non- 
rotating early-type galaxies the ability to retain massive halo of hot. 
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X-ray emitting gas corresponds to a distinct low-ionisation charac- 
ter of tiie emission-line regions and to a remarkable spatial coin- 
cidence between optical and X-ray emitting features. These find- 
ings testify to an interaction between the hot and warm phases of 
the ISM that provides an additional excitation mechanism besides 
photoionisation by old stellar sources. 

The prominent role of old stellar sources over other ionising 
mechanisms in powering the nebular emission of early-type galax- 
ies, taken together with the extended character that ionised gas al- 
ways displays in our sample, motivated us to explore the nature of 
the ionised-gas emission that is observed in early-type galaxies at 
larger distances through large physical apertures, as in the case of 
the SDSS spectroscopic data. We have exploited the integral-field 
nature of our data to integrate the SAURON spectra of our sample 
galaxies over increasingly wider circular apertures as to mimic the 
spectra of increasingly distant objects observed through the fixed 
3" -wide SDSS aperture. We measured the values for the [O III] 
equivalent width and the [O IIl]/H/3 ratio for the total emission 
in such apertures and compared our results with the emission-line 
measurements f or the morpholog i cally- selected sample of early- 
type galaxies of ISchawinski et aU ( |2007|) . Even accounting for the 
incomplete character of the SAURON sample, our data suggest that 
in very few, if any, of the SDSS galaxies which display only mod- 
est values for the equivalent width of the [O III] line (less than 
~2.4A) and LINER-like values for the [O IIl]/H/3 ratio, the neb- 
ular emission is truly powered by AGN activity. Only the most in- 
tense, and rarer, manifestations of LINER nuclear activity can be 
detected against both the stellar background encompassed by the 
SDSS aperture and the diffuse emission of early-type galaxies. The 
latter, when detected through SDSS spectra, is usually powered by 
pAGB stars associated either to old stellar populations or to a recent 
star-formation episode. 
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